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Abstract 

In a theory with broken supersymmetry, gaugino coupHngs renormaHze differently from gauge 
coupHngs, as do higgsino couplings from Higgs boson couplings. As a result, we expect the gauge 
(Higgs boson) couplings and the corresponding gaugino (higgsino) couplings to evolve to different 
values under renormalization group evolution. We re-examine the renormalization group equations 
(RGEs) for these couplings in the Minimal Supersymmetric Standard Model (MSSM). To include 
threshold effects, we calculate the /^-functions using a sequence of (non-supersymmetric) effective 
theories with heavy particles decoupled at the scale of their mass. We find that the difference 
between the SM couplings and their SUSY cousins that is ignored in the literature may be larger 
than two-loop effects which are included, and further that renormalization group evolution induces 
a non-trivial flavour structure in gaugino interactions. We present here the coupled set of RGEs 
for these dimensionless gauge and "Yukawa"-type couplings. The RGEs for the dimensionful SSB 
parameters of the MSSM will be presented in a companion paper. 
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I. INTRODUCTION 



Renormalization group equations (RGEs) have played an important role in extracting 
phenomenological predictions of theories valid at very high energy scales. One of the best 
known examples of this is the prediction of the weak mixing angle from the simplest SUSY 
SU{5) grand unified theories (GUTs): the equality of the gauge couplings at Q = Mgut, 
which implies sin^ OwiMQu^^) = 3/8, leads to the measured value of sin^ Oyy only when these 
couplings are evolved to the scale Q ~ 100 GeV relevant to experiments Since at least 
the 1980's, RGEs have played a role in the analysis of the implications of widely different 
models of new physics, including supersymmetric models, the subject of this paper 

Supersymmetry (SUSY) is a novel symmetry that links bosons and fermions, provid- 



ing a level of synthesis never before attained [3|, 14, Is, B]. Since supersymmetry implies a 
super-partner for every chirality state of the Standard Model (SM) with the same internal 
quantum numbers, it implies an approximate doubling of the degrees of freedom. Supersym- 
metry, even if softly broken, fixes the dimensionless couplings that determine the tree-level 
interactions of the superpartners with SM particles and with one another [i.e. interactions 
via mass dimension four), so that in this sense the theory is completely predictive. The 
interest in SUSY phenomenology was sparked once it was understood that the SUSY non- 
renormalization theorem implied the stability of the hierarchy between the weak and GUT 
or Planck scales to radiative corrections [3]. Since SUSY is clearly not manifested in the 
observed spectrum of elementary particles, SUSY must be (explicitly or spontaneously) bro- 
ken. Moreover, the mass scale of SM superpartners - at least of those that couple sizably to 
the Higgs sector - must be smaller than 0{1) TeV, in order for SUSY breaking not to desta- 
bilize the hard- won hierarchy between the weak and GUT scale just mentioned. In view of 
our complete ignorance of the mechanism of SUSY breaking, it is expedient to parametrize 
SUSY breaking by introducing soft-supersymmetry-breaking (SSB) operators with dimen- 
sion < 3 that respect Lorentz, gauge and any other symmetries into the effective Lagrangian 
of the theory. These include gaugino mass parameters, the SSB Higgs mass parameter b of- 
ten written as Bu, sfermion mass matrices, and trilinear scalar coupling parameters, usually 
denoted by Aijk |8|. 

We emphasize that the proliferation of SSB parameters is a result of our ignorance of how 
SUSY breaking is communicated to MSSM super-partners. Until this is understood, we are 
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led to resort to a variety of models [a, lid, I |ll|,ll^. Typically these models, which are based on 



differing assumptions about high scale physics, lead to different ansatze for the pattern of SSB 
parameters, in a Lagrangian renormalized at this high scale Qhigh- This Lagrangian cannot 
directly be used to perturbatively extract phenomenological predictions at the E ~ 100 GeV 
scale of our experiments, because large logarithms log (QhighZ-E') would invalidate (fixed 
order) perturbation theory. Instead, RGEs are used to evolve all the Lagrangian parameters 
';o the scale E, and this resulting "low energy" Lagrangian is used to evaluate its implications 

Since supersymmetric models automatically contain many scalars with the same gauge 
quantum numbers, there are potential new sources of flavour violation [?, [3]. Within the 
parametrization of the SSB terms described above, these are encapsulated in the sfermion 
mass matrices and in the A-parameters (which may be regarded as matrices in the flavour 
space). Our ultimate goal is to develop a program that will enable a general study of flavour 
violation, for arbitrary values of SSB parameters at Q = Qugh- For the present, we confine 
ourselves to quark flavour violation.^ 

Toward this end, we need to know the RGEs for all dimensionless and dimensionful 
parameters. These are already known, both in the SM and in the MSSM Q, Q, H Q 
at the two-loop level. Working to two loop accuracy means that we must include threshold 
effects in the one-loop RGEs, since these can numerically be as large as the two loop terms, 
and perhaps even larger if the sparticle spectrum is very split. It is in this connection 



that we extend the existing literature 19|, l20|| as described below. As in these studies, we 



implement SUSY and Higgs particle thresholds as step functions in our evaluation of the 
/3-functions. We consider first the evolution of the dimensionless couplings. We assume 
that for any particle with mass between the weak scale and the scale Qhigh introduced 
above, the effective theory used to calculate the /3-function includes this particle if Q > Mj, 
but not otherwise. Thus, at a high scale Q above the masses of all sparticles and Higgs 
bosons, the /5-functions are those of the MSSM. As we come down in scale, the SUSY 
particles and the heavy Higgs bosons decouple one- by-one until we are left with just the 
SM particles, and the /^-functions for the dimensionless couplings of the SM. Although 



^ The analysis can readily be extended to lepton flavour violation. Although lepton flavour violation may 
be more striking and experimentally easier to access than quark flavour violation, it may depend on yet 
other unknowns such as the structure of any singlet neutrino mass matrix. 
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this sounds quite straightforward in principle, there are several of issues that need to be 
confronted to implement this program. 

• The step-function decoupling of particles at the scale of their mass makes sense for 
couplings of fields only in their (at least, approximate) mass basis. For instance, in 
the MSSM Higgs sector, it is quite clear where the spin zero mass eigenstates h or 
H (or A and H^) decouple, but the decoupling point for the hu or hd fields (or of 
the higgsinos) is ambiguous, because the mass eigenstates are combinations of the hu 
and h(i type fields (or the corresponding higgsinos). Moreover, since the Higgs bosons 
and higgsinos acquire mass from different origins, these are not "diagonalized" by a 
common rotation, so that in order to implement the decoupling we must in principle 
work with Higgs boson/higgsino bases that are not SUSY transforms of one another. 
See, however. Sec. IIII Al for a simplification. 

• Since we decouple the superpartners but not the lighter SM particles, the effective 
theory below the scale of the heaviest sparticle "knows" about SUSY breaking, and 
necessarily includes new dimensionless couplings. For instance, while in the SUSY 
limit, the coupling of the gauge boson to quarks is exactly the same for the corre- 
sponding gaugino- quark- squark coupling, below the scale where we begin to decouple 
heavy particles, these couplings no longer run together. For every gauge coupling gfj, 
the low energy theory may include additional distinct gaugino-particle-sparticle cou- 
plings gf (labelled by the scalar $) that evolve differently from Qi. For instance, if 
ntq <C Trig, below Q = rrig the bino coupling g'^ evolves differently from g'. Moreover, 
g"^ and g'^ (as well as other couplings) also do not evolve in the same way. Indeed, 
as we will see later, these couplings also develop flavour-violating components and so 
become matrices in the flavour-space, which we will denote by gf , once again labelled 
by the scalar $. 

• The proliferation of these new fermion-fermion-scalar couplings also extends to the 
Higgs sector. For instance, if the heavier Higgs bosons H^, A and H all decouple at 
the scale Q ^ tuh that is hierarchically larger than and squarks are light, the 
effective theory in the range \fi\ < Q < niH includes both higgsino doublets /i„ and hd 
coupling to various flavours of quarks and squarks, but with quark (and squark) pairs 
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coupling only to the light SM-like Higgs boson h. Thus, we must also have independent 
higgsino coupling matrices that we denote by f^^e, once again, labelled by the scalar. 

We see that even for the dimensionless couplings, there are new complications which (to 
our knowledge) have not been included in previous studies, but must be taken into account 
for a proper implementation of threshold effects into the (one-loop) RGEs. Motivated by 
the fact that our goal is to study flavour violation in sparticle decays in a general way, 
we derive the RGEs for all the relevant dimensionless couplings of the MSSM, including 
flavour-violating superpotential interactions, in this paper. The number of RGEs for the 
dimensionless couplings now expands from those for the three gauge couplings and three 
Yukawa coupling matrices to those for the three gauge couplings, nineteen scalar-fermion- 
fermion coupling matrices and four Higgs-higgsino-gaugino couplings. We should mention 
that, in addition, there are numerous quartic scalar couplings arising from the D-terms or 
the F-terms in the scalar potential. Although these are given by the "squares" of gauge or 
superpotential Yukawa couplings in the SUSY limit, they will also renormalize differently 
below the SUSY thresholds. We do not discuss the evolution of these couplings which are 
less important phenomenologically. Fortunately, these couplings do not enter the evolution 
of the dimensionless couplings that we do consider in this paper. We derive the RGEs for the 
dimensionful parameters, the scalar masses (matrices in the case of squarks), the trilinear 
coupling matrices and the SSB Bfi parameter, which have their own complications, in a 



companion paper 



2l| 



Since, as we mentioned above, the threshold effects necessarily involve SUSY breaking, 
we derive our RGEs from the RGEs in a general {i.e. not necessarily supersymmetric) field 
theory that have been worked out in the literature. We use the seminal papers by Machacek 
and Vaughn jl^ to derive the RGEs for the dimensionless couplings, while we use the work by 



Luo, Wang and Xiao 22| to derive those for the dimensionful couplings of the MSSM. These 
studies both use the two component spinor formalism. In contrast, most phenomenologists 
are much more familiar with the four-component formalism. To facilitate the use of the 
general RGEs, we re-cast the formulae for the RGEs for dimensionless couplings into four- 
component notation in Sec. [TTl In Sec. IIIII we discuss the transition to the "mass basis" 
necessary for the implementation of the MSSM thresholds. The actual derivation of the 
RGEs is carried out in Sec. IIVI In Sec. |V] we show numerical examples of the solutions to 
some of the RGEs and discuss the impact of the new effects we have included in our analysis. 



We summarize in Sec. I VII The set of RGEs for the dimensionless couphngs is hsted in the 
Appendix. 



II. FORMALISM 



The one- loop RGEs for gauge couphngs are weU known and take the form 23 1 



-9 



S{Rs 



fermions 

where C{G) is the quadratic Casimir for the adjoint representation of the associated Lie 
algebra, and S{Rp) and S{Rs), respectively, are the Dynkin indices for the representations 
Rp, Rs under which the fermions and (complex) scalars transform. For the Lie algebra of 
SU{N), C{G) = N, while S{R) = 1/2 for the fundamental A^- dimensional representation, 
and S{R) = N for the adjoint representation. For the f/(l)y gauge coupling g', C{G) = 
while S{R) = {Y/2y. Heavy scalars and fermions are decoupled simply by excluding them 
from the sums on the right-hand-side of ([T]). 

The derivation of the RGEs for "Yukawa-type" couplings of scalars with fermions is more 
involved. As discussed in Sec. [Hour strategy for deriving the RGEs is to first cast the general 
two- component results in Ref. [l^ 



and Ref. 



22 



into four-component notation, match the 



parameters with those of the MSSM Lagrangian density, and write down the RGEs. Toward 
this end, we begin with the Lagrangian density written in terms of two- component spinor 
fields ^/'p and real scalar fields (pa, as 



C 



(2) 



ri'fipAi^ A 



1 

2 L 



h.c. 



^ab(Pa(pb 



-Y^,^/;JC^g0a + h.c. 



4 

1 

2! 

— h.abc4>a4>b4>c 



(2) 



The anti-symmetric matrix = i(T2, which is included to make the spinor bilinears Lorentz 
invariant, also ensures that the Yukawa coupling matrices Y"^ (labelled by the scalar field that 
couples to the spinors) are symmetric in the fermion field type indices p, q. The spinor mass 
matrices mj are likewise symmetric, but not necessarily Hermitian. The coefficients rriab, 
habc and Xabcd are real and completely symmetric under permutations of all their indices.^ 



^ The factor i in front of the Yukawa couphngs, which is a reflection of the symmetry of the Yukawa 
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We now turn to the corresponding Lagrangian density in the four- component notation 
with Dirac and Majorana spinor fields and complex scalar fields, usually used by phe- 
nomenologists. The familiar four-component Dirac spinor has unrelated left- and right- 
chiral components and is, therefore, made up by two independent two-component spinors ipL 
and ipR according to = {ipi, — C'^/j) (where we really mean to write this equation as one 
for the column matrix \1/d)- In contrast, the left- and right- chiral components of a Majorana 
spinor which, by definition, satisfies = C^\j are related. Thus, for example, a Majorana 
spinor may be written in terms of a single two-component spinor as = {"^Li ~C'^l) 
alternatively, as The general form of the Lagrangian density that we work with 

in the context of the MSSM with a conserved i?-parity then takes the form, 



Mk 



+ 



2! 



Bab^a^b + h.C. 



(U^)^., ^>o,PL^Dk<^a + (U^) ^D^PL^Dk^i 



jk 



+ 1 (X^)^.^ ^MjPL^Mk<^a + I ^MjPL'^Mk^i + h.C. 



+ 



2! 



^lliabc^b^^ + h.C. 



1 1 

2! 2! 



AabcdH'^l'^c'^d 



3! 



Kbcd'^i'^b^c^d + h.C. 



(3) 



In ([3]), Yukawa couplings of Dirac fields to scalars are denoted by the matrices and U^, 
while couplings of scalars to one Majorana and one Dirac field are denoted by the matrices 
Vq and Wq. These matrices have no particular symmetry (or hermiticity) properties under 
interchange of the fermion field indices j and k. These indices label the fermion field type 
(quark, lepton, gaugino, higgsino) and also carry information of fiavour and other quantum 
numbers {e.g. weak isospin and colour). On the other hand, the matrices and X^ that 
couple scalars to two Majorana fields are symmetric under j ^ k because of the symmetry 
properties [sl of the Majorana spinor bilinears that appear in ([3]). The scalar mass squared 
matrix is Hermitian in the scalar field indices, a, b, while the matrix B is symmetric. 
The Majorana fermion mass matrices nix and are both Hermitian and symmetric. The 

coupling matrices, was included in Ref. 22 1 but not in Ref. [l^ . 
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matrix mv is the familiar mass matrix, while mV is the coefficient of the CP- violatiner 

n 

Majorana fermion bilinear p|, whose effects in the two-component spinor language show up 
as phases in the entries of the fermion mass matrix mj in ([2]). The trilinear and quartic 
scalar couplings A, A' and H are symmetric under interchanges a ^ b and/or c ^ d for A, 
under interchanges of b, c, d for A', and finally under b ^ c for H. 

Before proceeding further, we should point out that the Lagrangian density in ([3]) is not 
the most general one that we can write. For instance, we have not included mass terms for 
Dirac fermions. Also, in writing the Yukawa interactions of fermions, we have not included 
terms with the operators 

"^a/Pl^d*, and ^/jPl^a/$^ . 

Dirac fermion masses appear in the MSSM only upon the spontaneous breakdown of elec- 
troweak symmetry, while the omitted Yukawa terms do not appear in the MSSM with 
i?-parity conservation because Dirac fermions then carry a baryon or lepton number which 
is assumed to be conserved by these renormalizable operators. Likewise, we have also written 
only the subset of all scalar interactions that we need for our analysis of the MSSM. 

Upon substitution of the four- component spinors in terms of their two-component cousins, 
and decomposing the complex field 



V2 

into its real and imaginary pieces, both of which are real fields, we can re-cast (jS]) into the 
general form of (I2l), involving only two-component spinors and real scalar fields. This then 
allows us to translate the parameters in into the corresponding parameters in ([2]). We 
can now use the general results for the RGEs for the dimensionless [l^ and dimensionful 



22l | parameters of a general quantum field theory to obtain the RGEs for the dimensionless 
Yukawa coupling matrices as well as for the mass and trilinear coupling parameters that 
appear in (I3l). We then use these to obtain the RGEs for the parameters of the MSSM, 
including the effects of decouplings of various sparticles. 

Since, our focus in this paper is on the RGEs for the dimensionless parameters contained 
in the matrices U^, U^, V^, Wq, and X^, we begin by re-casting the corresponding 
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terms into two spinor component notation to find, 

( (uf + uf 



(4) 



3 - 



\ 



( 











u 



u 



2T 



where 





V ^ 



ai 



(4) 



+ h.c. 



with the Dirac spinor = ("^l, — C"^/?) and the Majorana spinor \1/a/ = {ipM, —Ci^li)- 
Notice that the two-component spinor fields are coupled to real scalar fields, and 
via symmetric (but not necessarily Hermitian) matrices as expected. We can now use the 



one-loop RGE for the Yukawa matrices Y that appear in ([2]) 



14 



22 



(5) 



(47r)' /3?-li„ioop =1 [YUF)^" + Y'^Y^iF)] + 2Y^Y"tY^ 

+ Y^Tr {| (Y^tY'^ + Y""^Y^)} -3/{C2(F),Y"} , 

to obtain the one-loop /? functions for the various dimensionless couplings in ([3]). Here, 
Y2(-F) = Y^'I'Y^ and C2{F) = t^^t"^, where are the group generators in the reducible 
representation that includes all the fermion fields. Clearly, the matrices Y", Y2(-F) and 
C2{F) all have the same dimensionality, determined by the total number of two-component 
fermion fields in the system. We also draw the reader's attention to the fact that in the 
four-component notation the sum over a, 6, ... runs over the different complex scalar fields, 
whereas in the two-component notation, we must not only sum over the various field types, 
but also separately over the real and imaginary parts of the same complex field. In other 
words, we have two distinct Yukawa matrices Y", one when a refers to and the other 
when a refers to ^ai- The trace that appears above is a sum over the fermion types. 



Eq. ([5]) is slightly modified from that in Ref. Q ■ We have written the (F) instead of Y2 {F) in the 
first term and symmetrized the trace with respect to a and b. The second modification also appears in 
Ref. 



22l | , while the first one preserves the symmetry of the Yukawa coupling matrix. 
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The /3-functions for the Yukawa couphng matrices in ([3j) can then be readily found to be, 



(4vr)'/3f/i|,. 



loop 



(4vr)'/3f/.|,. 



'b^b 



+ U^Tr { (Ufui + U^U^) + i (X^txi + xfX^) } 
- 3/ KC^(F) + Cf (F)Ul] 

+ Iv, (vtU^ + Xf W,) + 2U^Wt W, + U.^Tr { WtW, + Vtv,} 



(6) 



loop 2 



+ U,^Tr { (V'M + Vl^Vl) + I (Xl^^l + X^tx^ 
+ U^Tr I [V'M + U^tu^) + i (x^tx^ + xitx^ 
- 3/ [U^C^(F) + Cf (F)U^] 

+ ^ (U^WJ + V,X,^t) W, + 2V,VtU^ + U^Tr { Wjw, + Vlv,} 



T2tTT2 



(7) 



Within the MSSM with i?-parity conservation, the last line of each of these equations van- 
ishes because and vanish for the a values for which and are non-zero. Con- 



10 



tinuing, 



(4^)'/3^ili-ioop=^[(w.Wt + V^V, 



x^xfUi 



xi ( wtwT + v);v, + x^^x,^ + x^^x? 



+ 2 

+ x 



2(Xrxi 



x^tx^ 



x,^ 



Tr I (Uf + Vl^Vl) + I (Xf Xi + Xf X,^) } 



- 3/ [X^C^(F) + C^{F)Xl] 



+ 



+ 2 [W.Wtxi + XlWlWl] + X^Tr {wtw, + ViV,} 

(47r)^ /3x.Hi_ioop =^ [(W.WJ + V^V* + X^X^^ + x^xf ) X^ 
+X2 (w*w^ + VjV, + x^tx^ + Xf X^)" 
+ 2 [w.u^t V, + v^ur + x^x^tx^ + x^x^tx^] 



+ X^Tr 



2(X^^X^ 



- 3/ [X^C^(F) + C^{F)Xl] 



+ 



(v^u^* + x^w*) wj + (vlW, + w^x 

+ 2 [VlVlXl + XgVt V,] + X^Tr { W+W^ + Vt V, 



(8) 



(9) 



The last two lines of each of ([8]) and (Q vanish in the i?-parity conserving MSSM, again 



11 



because W° and are zero. Finally, 

(4vr)^ /3yJi_,oop =\ [(U^U^^ + U^Uf + V.Vj) V, 

+v, (w* + v;v, + X^txi + Xf X^^ 
+ 2 [Ui WtX^ + U^Wtxi + V,W:W^] 
+ V,Tr { WtW, + VjV,} - ?>g^ [V,C^(F) + C^(F)V,] 



+ 



wt + v,x^t + x^) + (u^ + u^) u,^t y 



wlx^ 



+ 2 (U^t + Uf ) V, + V, (X^t + Xf ) X^] 



+ VbTr 



u^^ (ui + u^) + (u^t + ) u: 
1 



+ 2 



(x^, + xa + (x^t + xf ) x: 



(10) 



(4vr)'/3H.Ji_ 



loop 



1 

2 LV 



f w,w; + v^v* + x^x.^t + x^xf ] w„ 



2 [V^V:W, + X^VtU,^ + X^VtU,^] 

W,Tr I W;W, + VtVfe} - 3(7^ [W,C^(F) + C^(F)W,] 

1 



+ 2 



w,uf + x^vn (u^ + u^) + (x^ + x^) v^u^ + x.^tw. 



+ 2 [W, (U^t + ) + X,^ (X^t + x^t) w,] 



W.Tr 



Uf (U^ + U^) + (U^t + uf) 



x^^ (x^ + x^) + (xi^ + xf)x^ 



2 L 

(11) 

The last four lines of each of ffTOj) and f|TT]) do not contribute in the -R-parity conserving 
MSSM because each of U^, U^, X^ and X^ vanish for the values of a for which Vq and Wa 
do not. 

The quantities C^{F) and C^(F) that appear in Eq. (l6l)-( fTTl) are the quadratic Casimirs 
C2{F) separated into contributions from left-handed, right-handed and Majorana fermions, 
so that 

C,{F) = diag (C^(F), C^(F), Cf (F)) . 
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The traces that appear in (l6|)- (fTT]) denote a sum over fermion types, which sometimes, but 
not necessarily, are the trace over fermion flavours when we refer to the MSSM. Note also 
that in the context of the MSSM, the index a may itself include a flavour (or even colour) 
index, for instance when = c^. 

The reader may legitimately wonder what we have gained by writing the RGE which took 
the much simpler- looking form ([5]) in the much more cumbersome-looking and definitely 
much longer forms given by Eq. fl6l)- f|TTl) . The reasons for doing so are two-fold. First, many 
authors are more familiar with the MSSM couplings written in four-component notation, so 
that the longer four-component equations can be more directly used. Second, and perhaps 
the more important reason is that much of the work that would he needed to obtain the 
MSSM RGEs, even starting with the couplings in two-component notation, has already been 
done when the RGEs are written in as in W\)- ni\) . The many terms that appear here are 
indeed present in the MSSM RGEs (as we will see) and starting with the compact notation 
of (IS]) only means that part of the work has to be repeated each time we use it to get the 
RGE that we need. 

Before closing this section, we remark that the RGEs for the dimensionless gauge cou- 
plings and those for the dimensionless "Yukawa-like" couplings form a closed set, uncoupled 
to the RGEs for the dimensionful parameters. Dimensionless quartic couplings of scalars 
also do not enter these RGEs. Since we have made no assumption about supersymmetry, 
this is true even with threshold effects included.^ In the next section, we turn to a dis- 
cussion of how threshold effects are to be included in the MSSM RGEs. Of course, since 
we want to work to two-loop accuracy, it suffices to include the threshold effects only in 
the one-loop terms, which is the reason for our focus on the one-loop RGEs of a general 



l4 and 



which is 



quantum field theory. We mention here that the authors of Ref. 
our starting point, use the MS renormalization scheme rather than the DR scheme better 
suited for the analysis of radiative corrections in supersymmetric theories. We can, however, 
use their general expressions (only for the one-loop part) to obtain our RGEs since there is 
no scheme- dependence at this level. We can then augment these with the two-loop terms 



(without threshold corrections) of the MSSM RGEs 



15l | to obtain two-loop accuracy in the 



DR scheme. 



Of course, in this case, information about the SUSY spectrum enters via the focation of the thresholds. 
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III. PARTICLE DECOUPLING 



We are now equipped to derive the RGEs for the dimensionless couphngs in any quantum 
field theory with an arbitrary set of Dirac and Majorana spinor fields together with spin zero 
fields (with non-trivial gauge quantum numbers) whose Lagrangian can be written in the 
form ([3]). This is, as we have already noted, not the most general form for the Lagrangian 
density since some field operators, such as mass terms for Dirac fermions or dimension 
four operators of the form '^mPl'^d^ and dPl'^ that do not occur in the MSSM 
with i?-parity conservation imposed have been omitted. Since ([3]) makes no reference to 
supersymmetry, we can continue to use it even when SUSY is broken via the inclusion of SSB 
terms. The procedure for obtaining the one-loop RGEs is now straightforward. We simply 
write down the Lagrangian density for the MSSM, including the SSB terms consistent with 
Lorentz and MSSM gauge symmetries and the assumed conservation of i?-parity, compare 
the dimensionless couplings with those in ([3]) to construct the matrices U^, U^,- ■ ■ ,Wa, and 
then use (l6l)- (fTTl) to get the required RGEs. This program has two complications that we 
must take care of in its implementation. Although we have referred to these in Sec. [H it 
does not hurt to repeat them here. 

• For Q values above all SUSY thresholds, the effective theory is supersymmetric, and 
couplings of SM particles are related by SUSY to corresponding sparticle couplings, 
so that the particle and sparticle couplings evolve the same way. This is, however, not 
necessarily the case below the SUSY thresholds. Then, while the /3-function for say 
the hypercharge gauge coupling still depends only on the hypercharge gauge coupling 
(because of Ward identities, its one-loop evolution is governed by just the "vacuum po- 
larization corrections"), the /^-functions for the corresponding "Yukawa- like" couplings 
g'-^ of the bino to the fermion-sfermion pair will, in general, depend on gauge as well 
as all "Yukawa" couplings, including those from super-potential interactions. There 
is no principle that precludes this once supersymmetry is broken. However, above 
all SUSY thresholds, this P function must reduce to that for the gauge coupling, so 
that all dependence on super-potential Yukawa couplings must cancel. This provides 
a non-trivial check on our formulae. Thus our system must be extended to include 
the couplings gf as well as the couplings f * the additional "Yukawa-like" couplings of 
gauginos and higgsinos, respectively, which as we have noted become matrices in the 
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flavour-space. 

• To implement SUSY and Higgs particle thresholds, we need to have an idea of the 
spectrum as well as the couplings in the corresponding mass basis. Of course, the exact 
spectrum is only obtained after extracting the weak scale values of all SSB parameters, 
and then diagonalizing the various mass matrices. Since the positions of the thresholds 
only enter logarithmically in the solutions to the RGEs for the various gauge and 
"Yukawa-like" couplings, it suffices to include these in a reasonable approximation. 
Obviously, this entails rewriting the interactions in terms of the approximate mass 
basis for the fields. A discussion of how we do so forms the subject of the remainder 
of this section. 

A. The Mass basis and MSSM thresholds 

In addition to the sparticles of the SM, the MSSM includes squarks and sleptons, gluinos, 
charginos and neutralinos, and the additional spin zero particles A, H and in its Higgs 
sector. At a scale Q larger than the masses of all these particles, the RGEs are given by 
those of the MSSM, and as we reduce Q, these morph into those of effective theories that 
interpolate between the MSSM and the SM, and ultimately at the scale Q ~ m^, reduce to 
those of the SM, assuming that all sparticles and additional Higgs bosons are significantly 
heavier than m^. This will serve as another check of our formulae. 

Assuming for the moment that /l and fn are also approximate sfermion mass eigenstates, 
i.e. that mixing among the sfermions is modest, decoupling the sfermions is straightforward. 
Simply exclude any sfermion field (and all its couplings) in the effective theory at a scale 
Q below its mass, and evaluate the RGEs using this truncated theory. The same is true 
for gluinos which, being colour octets, cannot mix with any other fermion since SU{?))c is 
unbroken. The decoupling of charginos, neutralinos and the additional Higgs bosons is more 
complicated because these are frequently strongly mixed. 

Our procedure for decoupling these fields is guided by the observation that threshold 
effects that we are trying to evaluate are significant only if there are several well separated 
mass scales Mi in the spectrum. If this is not the case, then we can decouple all the additional 
particles at a common scale Msusy, since any error would only be of O (^-^ In ^^^^ and 
would ^ as Mj ^ MgusY- 
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In the Higgs boson sector, threshold effects are, therefore, important only if the mass 
scale of the additional bosons ^ rrih, in which case rriA — rnn — rnH± ^ rrih < rrit. We can 
then decouple all the additional spin- zero bosons at a common scale that we take to be rriH, 
below which h plays the role of the SM Higgs boson. To implement the decoupling, we must 
rewrite the Lagrangian couplings in terms of the (approximate) mass eigenstates contained 
in the SU{2) doublets^ 



(12) 










) ^\ht 
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\h': 



(13) 



where the electrically neutral, complex fields h and Ti are given by. 



h 



H 



h+iG° 

-H+iA 
V2 



and s = sin/9 and c = cos/3. Above Q = tuh, the theory includes two Higgs doublets which 
we may take to be hu and hd, or any two orthogonal combinations of these. However, at 
Q = niH, it is specifically the doublet (IT^ that decouples, and the effective theory includes 
just one spin-zero doublet which couples to quarks and leptons via Yukawa coupling matrices 
Xu = fuSin/3, Xd = id cos (3 and Ae = feCos/5, with these relations imposed at Q = mn- 
For Q < tuh, we will be interested in the evolution of the Aj's which are the couplings 
of the doublet f|T2l) that remains in the theory. We stress again that supersymmetry is 
broken, and the decoupling of the second doublet makes no statement about the decoupling 
of the corresponding higgsinos. The effective theory below Q = rriH may or may not contain 
higgsinos, depending on, as we will see next, the value of /i. 

In the chargino-neutralino sector, threshold corrections are small if \fi\ and the SSB 
gaugino masses all have a magnitude ~ Mz, since their effects will then persist down to the 
weak scale. Threshold effects are significant only when at least one of these is hierarchically 



^ We use the notation of Ref. [3|] where the left-chiral up-type Higgs superfield iJ„ transforms as a 2 of 
SU (2), while the left-chiral down- type Higgs superfield transforms as a 2*. The doublets that make up the 
linear combinations in and both transform as a 2 of 5[/ (2) and have a positive weak hypercharge. 
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larger than the weak scale, in which case gaugino-higgsino mixing effects are small, and the 
gaugino-higgsino mass eigenstates can be approximated by the bino, the wino-triplet and 
the two Majorana higgsino doublet fields,^ 

K2 = V^'^ . (14) 

To implement decoupling, we must rewrite interactions of gauginos and higgsinos in the basis 
consisting of binos, winos and these two higgsino doublets. Above the scale of the higgsinos, 
we could have worked with any two orthogonal linear combinations of the higgsino doublets, 
whereas below Q = \fi\, both higgsinos decouple from the effective theory. We note here that 
we have to perform different rotations in the spin-zero and spin-| Higgs sectors to go to the 
mass basis. This appears to be different^ from what has been done in the literature 19|, |20 |. 
There is one additional complication that enters when we rewrite the Lagrangian density in 
the (approximate) "mass basis" . Some of the fermion fields have negative eigenvalues, and we 
are led to make field re-definitions [sl ^ {i'y^y'''^k, where 6^ = 0(1) if the corresponding 
mass eigenvalue is positive (negative), on Majorana spinors to get their masses in canonical 
form. While this introduces additional 75 's and phases in the interactions, we have checked 
the RGEs are independent of these, as may be expected since these RGEs do not depend 
on the masses. 

Finally, let us return to mixing effects in the sfermion sector. We have ignored these in our 
analysis. For the squarks of the first two generations, flavour-physics constraints typically 
restrict the size of inter-generation mixing effects (in the basis where the corresponding 
quark Yukawa couplings are diagonal). There is, however, no principle that dictates intra- 
generation mixing to be small, even though this is usually assumed to be the case for first and 
second generation fermions in many models, while third generation sfermions are expected to 
mix via their Yukawa couplings*. This said, barring large, accidental cancellations we would 
expect that the physical sfermion masses are approximately of the same order, in which case 



^ and tphd are the fermionic components of the SU{2) doublet left-chiral superfields Hu and Hd- 

Following Ref. , we denote the charged and neutral components of tph^ by (^/'^+ , tphl ) and of iph^ by 
{ijjf^- ,tph°)- is a Majorana spinor whose left-chiral component destroys positively charged higgsinos, 
while its right-chiral component destroys negatively charged higgsinos, and similarly for -0,^- . 

_ d 

' We mention though that because we choose to decouple both higgsinos at the common scale |^|, this 
rotation of the higgsino fields is irrelevant for practical purposes. 

^ Substantial intra-gcncrational mixing may occur via the trilincar sfcrmion-fcrmion-Higgs boson SSB pa- 
rameter which is, after all, independent of the corresponding super-potential Yukawa coupling. 
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threshold effects from intra-generation mixing are not important, or the diagonal SSB terms 
are themselves hierarchical so that the unmixed sfermions fi and fn are approximate mass 
eigenstates, and can be decoupled at their physical masses. 



IV. APPLICATION TO THE MSSM 



A. Interactions 



We now use the RGEs for the general quantum field theory that we obtained in Sec. [ITl 



to derive the RGEs 
the notation of Ref. 



br the dimensionless couplings of the i?-parity conserving MSSM. In 
b|, the superpotential for the MSSM is given by, 



while the SSB terms may be written as, 

>CssB = - |M[fc(mJ)MML, + d\^{m^Q)kidLi + v}R^,{ml)kiUm + d}j^^^{my)kidm 
+ 4fc(mi)fc/^Li + e[k{ml)kieLi + 4fc(m|)fc/era 

+ {uRki^DkiULihl - u^^^^i&DkidLiK + dRki^d)kidLihd 

(16) 

- i [MiAoAo + MaApAp + Mgl^^^] 

% — — — 

- - [M{Ao75Ao + M'^Xp-f^Xp + M'^gA-^^QAl 

+ {h{Kh- + hlhl) + h.c.]. 

where the matrix indices i, j, ■ ■ ■ label the (s)particle flavours. The sfermion mass squared 
matrices itiq, m^, . . . are Hermitian, whereas the Yukawa coupling matrices fn,d,e that 
appear in the superpotential, as well as the trilinear SSB scalar coupling matrices a^^^^e have 
no hermiticity or symmetry property. Finally, the gaugino mass parameters Ma and are 
real, while the SSB Higgs mass parameter may, in general, be complex. Since the RGEs 
for the dimensionless couplings do not depend on any dimensionful parameters, we do not 
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really need £ssb for their derivation. The SSB parameters only enter the analysis in that 
they determine the positions of the various thresholds. 

In order to be able to compare our results (above all thresholds) with the standard 
results already in the literature 1^, 16, [itI. [isl we need the following conversion between 

isl: f = Y^, a = -h^, b = -B. Moreover the 



the notation used here and that in Ref. 

complex gaugino mass parameters (Ma) in Ref. [15| are equivalent to Ma — iM'^ in (fT6|) . 
It is clear that in order to obtain the RGEs, we need to extract the various matrices 

12 12 

U^' , X^' , V$ and W$, whose elements are the couplings between the various scalars, $, 
and fermions of the MSSM. Thus fermion-fermion-scalar interactions play an important role 
in our derivation. The interactions of quarks and squarks with Higgs bosons and higgsinos 
may be obtained from the superpotential, and are given by, 

+ -e,{Qlh-,PLU, + e,(fe)J,/i°Pie, + h.c] 

+ ^^-~e%(S:%PLU, + e,{i^)luL^PL^^- + ^^o^e%{il%PLe. + e,(fi^)J,eL.PL*/.o 
+h.c.] 

(17) 

We have written these interactions with a left chiral projector sandwiched between the four- 
component spinors to facilitate comparison with Eq. ([3]). Notice that the couplings of matter 
fermions to Higgs bosons are denoted the same way as the superpotential Yukawa matrices. 
In contrast, we have denoted the couplings of the fermion-sfermion system by a matrix with 
a twiddle on it, labelled by the scalar that enters the corresponding interaction. Above 
the scale of all particles and Higgs boson thresholds, the effective theory is the MSSM and 
these f* couplings are the same as the corresponding superpotential couplings f , but may 
be different below this scale. For instance, if < itlh-, the low energy theory just below 
Q = rriH contains two higgsinos, but just the one scalar doublet in (1121) . The effective theory 
below Q = rriH then contains only the couplings Xu,d,e along with the couplings f^^e; but 
no iu,d,e- As we will see the evolution of f*^g differs from the "would-have-been" evolution 
of iu,d,e for this range of Q. 

The remaining fermion-fermion-scalar interactions are those with the gauginos, and are 

19 



given by, 



where we have defined the following matrices: 







1 +h.c.| 







(18) 



I :g%-.A3 + i(g'%Ao (g%(Ai-zA,) ^^^^ 
(gQ),. (Ai + ZA2) -(g«),.A3 + i(g'«),.Ao 

^ I (g'^)..^3 - (g'"),.Ao (g^).. (Ai - 1X2) ^20) 
(g^),, (Ai + 1X2) -(g^),.A3 - (g'^),.Ao 

^'^"A3 + ^"^"Ao ^""(Ai-^As) \ 
^'^^ (Ai + zAs) -^^"Aa + ^'^-Ao / 

-^'"'^Aa - -g'^'^X^ ^^'<* (-Ai - zAz) \ 
G/.. = I • (22) 

^'^^(-Ai+zAa) ^^^Aa-^'^'^Ao / 

Here, Aq and Ap denote the (Majorana) hypercharge and 5'f/(2) gauginos. Just as the 

higgsino couplings in fll7p may differ from the corresponding Higgs boson couplings, below 

the various thresholds the effective theory is no longer supersymmetric, and the couplings 

of gauginos may deviate from the corresponding gauge boson couplings. Moreover, unlike 

gauge couplings, the renormalization of these fermion-fermion-scalar interactions depends 

also on the couplings f and f*, so that the g*'s can acquire a non-trivial flavour structure, 

as allowed for in (fTS!) above. Above the scale of the highest thresholds, however, the RGEs 

for these g*-type couplings must coincide with those for the gauge couplings (in the matrix 

sense) . 
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B. Gauge Coupling RGEs 



The one-loop RGEs for gauge couplings (including threshold effects) are well known. We 
will, however, write these for completeness. As explained in Sec. HI we will decouple particles 



using the "step approximation'' 
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20|. Toward this end, for any particle P, we define 



9p= 1 if Q > Mp, 
HQ <Mr . 

Throughout this paper, we will assume that Q > nit so that the SU (2) gauge symmetry is 
unbroken in the effective theory, and that there are three matter generations. Applying ([T]) 
for the MSSM particle content, and using 

3 



i=l 



leads to the familiar gauge coupling RGEs, 



4 + —Nn + —Ni,^ + —N-j + —Nf 
30 ^15 "^15 10 ^ 



5 ^ 



5 



22 1 1 1 

h 4 + -Nn + -A^f + - 

3 2 ^ 6 ^ 6 



3 ^ 6 ^ 



6 



26, 



11 



(23) 



(24) 



(25) 



Here, gi is the scaled hypercharge gauge coupling that unifies with the SU (2) and SU (3) 
couplings when the MSSM is embedded in a GUT: g'"^ = ^gf. We mention that although we 
have shown two distinct higgsino thresholds above, since we are working in the approximation 
that both higgsinos have the mass \fi\, we will from now on write 



Or 



C. RGEs for Yukawa and Yukawa-type Couplings 



Turning now to the RGEs for fermion-fermion-scalar couplings, we see that there is a 
large number of such "Yukawa-type" couplings in the Lagrangian. These are: the usual 
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Yukawa couplings to Higgs bosons, f^, f^, and fg; the couplings of higgsinos to the various 
fermion-sfermion pairs, f^, f^, ff', f^^, f^^, f|^; hypercharge gaugino couplings, g"^, g'"^, 
g'"«, g"^«, g''^«, the SU{2) gaugino couplings, g*^, g^, g^"", g^<^] and finally, the 

gluino couplings, g^, g"^, gf^. While the RGEs for the gauge couplings still form a closed 
set, those for the usual matter fermion Yukawa couplings do not. It is only the RGEs for 
all these fermion-fermion-scalar couplings that together with those for the gauge couplings 
form a closed system. These new couplings must coincide with the corresponding gauge or 
usual Yukawa coupling in the SUSY limit, i.e. above the scale of the highest threshold. This 
serves as the boundary condition for the RGEs for these new couplings. 

To implement decoupling of particles as described in Sec. IIIH we must construct the 
matrices U|,, U|, V$, W$, X|, and X|, for all $, in the (approximate) mass basis for the 
various fermion fields. We choose this basis to comprise of the Dirac fermions {uj, dj, z/j, Cj} 
together with the Majorana fermions {h\., h^, hf, hf, Aq, Ai, A2, A3, qa}- Here, ^ and h^^ 
are the neutral and charged components of the Majorana spinor h\^2 introduced in (fT4l) .^ 
The subscript % (which runs from 1 to 3) for the Dirac fermions is a fiavour index which can 
be suppressed if the Yukawa terms are written as matrices in fiavour space. This means that 
in the MSSM, and U| will be (4 x 4) blocks of (3 x 3) matrices when $ is one of the 
Higgs bosons in (fT2|) or (fT3|) . Similarly, since fiavour is carried by the sfermion scalar index, 
$ = /j, will be a (4 X 9) matrix where the number of rows can be further expanded to 
show each of the three matter fermion fiavours. Likewise, Wj. a (9 x 4) matrix where now 
the number of columns can be similarly expanded to exhibit these fiavours. Thus, when 
fully written out, Vj^ is a (12 x 9) matrix, while Wj. is a (9 x 12) matrix; see ( !29l) - (l32ll . 
Finally, and X| will both be (9 x 9) matrices. 

Within the MSSM with i?-parity conservation, the matrices U^'^ and X^'^ are non-zero 
only for $ = {h, 7^, G^, i/"*"}. They can be readily worked out by comparing (j3]) with the 
MSSM Yukawa interactions in ( fTTl) . It would take up too much space to display each of 
these matrices explicitly, so we will only present one example for each type of matrix. If 

^ As already mentioned, since we are working with a eommon threshold for higgsinos, we could have stayed 
in the original basis. We will see below that performing this rotation (by an arbitrary angle) gives us 
another check on our procedure. We also draw attention to the fact that we have combined the higgsinos 
of the 2 and 2* representations, so that although *re Majorana spinors, their chiral components do 
not annihilate a definite sign of the charge. 
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$ = h we have, 



/o \ 

cfj 



(26) 





V cf J / 

where we use a bold-type to indicate a (3 x 3) matrix in flavour space whose entries are 
all zero. The corresponding Uj!, matrix has its only non-vanishing entry in the (1, 1) block. 

2 
h 



The matrix X? takes the form 



X 







(4x4) 



ATf 



(27) 



(5x5) 



where 













gghu 


o\ 












—sgf^^ 












—isg^"" 








I 







isg^'^ 





0/ 



(28) 



The eight X|, and X| matrices can be obtained from flTSl) in conjunction with fl^Tl) and 
The reader will need to write these as symmetric matrices using the symmetry properties of 
Majorana spinor bilinears. Note that all the entries in the last row and last column of these 
matrices are zero because the Higgs boson has no coupling to gluinos.^° It is also clear from 
the MSSM Lagrangian that X^ contains only the g"^"^ and g'^'^ couplings, while X| contains 
only g''^'' and g^''. 

The V$ and W<j, matrices, on the other hand, are non-zero only for $ = 
{uL,dL,eL,h'L,UR,dR,eR}. From the interactions in flTTI) and (ITS]) we can see that for left- 
handed sfermions, Vj contains just f*-type couplings while contains only g*-type cou- 
plings. The opposite is the case for the right-handed sfermions. In the case where f = uu, 
we have 



V- 








(6x5) 



0(6x4) 0(6x5) 



(29) 
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We have suppressed the eolour index which would otherwise expand the last row/eolumn. 
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with 



V 



iiSZ (iSZ 0(3x1) 



(3x1) 



0(3x1) 0(3x1) (f^)^i (f^)^^ 



(30) 



As alluded to above, since the scalar labels on and carry sfermion flavour, we need 
to be careful when writing these if we choose not to expand out the matter fermion flavour. 
Here we denote this suppressed quark flavour index by the •, so that (f^)^j is really a (3 x 1) 
column matrix (since the index i is fixed by the squark flavour). Similarly 



W- = 



0(4x6) 0(4x6) 



Wu,. 



(31) 



(5x6) 



with 



i(g'^ 



0(1 



(1x3) 



0(1x3) 
0(1x3) 



2(%^)(g, 



s 



0(1x3) 
0(1x3) / 



(32) 



Notice that the positions of the suppressed fermion index and the (fixed) sfermion flavour 
index, i, are swapped between the matrices (|5Ui) and (15^ . Once again, we have suppressed 
the colour index on the quark (and squark) field. Also note that the entry from couplings 
to the gluino enter only via the matrix for ql, and only via the matrix for g^j, for 
all flavours of squarks. 

We are now ready to proceed with the derivation of the RGEs. Before we begin our real 
derivation, we perform a pedagogical exercise to facilitate comparison with earlier literature 
and also provide checks on our procedure. Since earlier studies did not make the distinction 
between couplings of quarks to gauge and Higgs bosons (or g^s and f 's) and the corresponding 
couplings (g's and f 's) of the quark-squark system to the gauginos and higgsinos, we first 
derive the RGE for (sf^) without separating out f 's and g's from f 's and g^s, respectively. It 
is also instructive to illustrate the additional complexity that results from keeping distinct 
thresholds for A, H and if^, and an arbitrary rotation angle for the higgsinos instead of 
the 45° rotation used in ([HI); i.e. 



*hO = (275)'HcVl/,o-s'^hO^ 



(^75)'^(s'vl>,o+c'vl/,o 



(33) 
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and 

where s' = sin (3' (set equal to (3 in Ref. [l^) so that the rotation is by an independent angle 
to the scalars. The factors (^75)^^'^ allow for the possibility that either higgsino could have 
had a negative mass eigenvalue, as noted below f[T^ . As we saw earlier, the RGE is not 
dependent on where or whether we include the (^75)^ factors. We stress again that we do 
not mean the equation that immediately follows to be the correct RGE, but write it only 
to facilitate comparison, and to illustrate some issues. Ignoring the separation of twiddle 
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terms on the couplings in the Lagrangian we find, 



(4vr) 



dt 



J + + 60+) + -{eH + eA + eH+\ 



kj 



kj 



+ s [s' {Oh - Ogo) + c2 {Oh - 6^)] {UXh 



+ s%),^ [3 (s^e, + c^^h) Tr{ftfJ + 



6h) Tr{3f]f, + ftf J 



3 2 Li! _ 3 ( 1^ + ^0 ) ^ 
5^^ 1 12 I 36 9 " ' ^ 



s 
2 L 



[s's + c'c)^e^o + (s'c - c'sfe 



en 



+ I {s'c - c's) {s's + c'c) {e~,o - e~,o) Oh 

s' {s's + c'c) — C {s'c — c's) 6}fi 



+ 



s's + c'c) 6*^0 + (s'c — c's) 6'^o 



4 "4''Q.- V2 . 
-2 (s'^s^ + c'^c^) + 2 (s'^c^ + c'^s^) 



+ 



[s'c — c's) (s's + c'c) ^6'^o — 6* 

-2sc {s" - c") (e,. - e,. 



1H 



s (s's + c'c) - c' (s'c - c's) ^ 



-2S ( S'^dr± + C'2^r± 



(35) 

A sum over repeated flavour indices k,l, . . . is imphed including over those with three re- 
peated indices, one of which is a 9q^,, e.g. OQi_{iuil)ik{^u)kj ^^^^ t'wo. Several comments 
are worth noting. 

• In writing the RGE as in (l35l) where we have ignored the difi^erences between the usual 
gauge/ Yukawa couplings and their twiddle counterparts, we have retained distinct 
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thresholds for each of the Higgs bosons as well as for the higgsinos. This is partly to 
facilitate comparison with Ref. and partly to indicate the origins of the various 
terms. 



We see that if we take the MSSM limit, where a 
we recover the usual, well known MSSM RGE 15 



1 the ^'s are set equal to unity. 
This is evident if 



16 
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the rotation angle that defines the doublets in ( fT2|) and ( |T3l) is treated as a scale 
independent parameter. In this case, the right-hand-side of (1351) is proportional to 
sin/5 which cancels out leaving the MSSM RGE for f^.^^ Notice also that in the 
MSSM limit, the dependence on the higgsino rotation angle, also disappears as it 
must. 

If, on the other hand, we set all the ^'s other than 9^, 6go and 6g+ to zero, the theory 
includes just the SM fermions, gauge bosons and the single Higgs doublet (fT2|) . that 
couples to SM quarks and leptons via the coupling matrices = sin (3, = fd cos j3 
and Ae = fe cos (3. It is then easy to see that we recover the relevant RGE for the 



Yukawa couplings in the SM 



241 ]. Indeed, all dependence on the f's in (15^ disappears 



in favour of dependence on the SM A's. 

Since SU{2)l remains a symmetry of the theory at low energy, we should expect that 
the couplings related to one another by SU(2)l should have the same RGE. We have 
checked that this is indeed the case, but only if we set a common value for all the heavy 
Higgs boson thresholds, and an independent common value for the higgsino thresholds. 
This should not be surprising since any splitting between the extra Higgs bosons, or 
between the higgsinos, is an SU{2)l breaking effect (and would entail introducing 
even more couplings into the low energy theory). We do not regard as reliable those 
terms in (1351) that come from SU{2) breaking effects because we are then using an 
inconsistent approximation. 



If instead we consider scale-dependent rotations, then the rotated fields themselves have an explicit scale 
dependence. The RGEs in Ref. [3, that were our starting point assume, of course, no such scale 
dependence of fields. However, we know that above all thresholds, the rotation (which is merely a field 
redefinition) cannot change the RGEs. In other words, there must be additional compensating terms in 
the RGEs which effectively allow us to take the sin/? out of the derivative on the left-hand-side, so the 
RGE once again reduces to the MSSM RGE. 



27 



Below the mass scale of the heavy Higgs bosons we may expect that since there is just 
one Higgs boson doublet in the theory, only the Xu4,e couplings enter. We can easily 
see that this is not the case if higgsinos are lighter than niH. There will still be two 
higgsino doublets in the low energy theory, and these couple via the matrices f^^e 
(notice the twiddle on the f). 

We have compared ( l35l) u sing a common threshold for the higgsinos with the corre- 
sponding RGE in Ref. 19|. We find agreement for all but the terms involving SU {2)^ 
and U{1)y gauge couplings. We find. 



^ ' dt 



-sit, 



CO 


fir / 




112 I 



1 n 4 

36 9 



+9l \ I - l%^w 



-6h — 36r,, 9to9 



to be contrasted with [1 



-s(fj [L2 /IZ + 3 _ / ]_ 4 1 \ 

V "Aj i^^yi 1 12 4 ''1 V36 9 ' 4 ^ 
f 9 9 3 / \ 



There is a similar point of disagreement in the RGEs for the down-type Yukawa 
RGEs.12 



• We found it difficult to perform a corresponding comparison with Ref. 20(] where there 
the RGEs appear to be written without doing any rotation of the Higgs fields, so that 
we could not abstract the relationship between the Higgs bosons /i. A, H and if^, and 
the thresholds that appear in their RGE. For the same reason, we could not see how 
to reduce this RGE to the SM RGE below the scale of all sparticles and heavy Higgs 
bosons. The MSSM limit is, however, correctly obtained. 

We reiterate that (135|) is not the correct RGE to use. If we derive this RGE keeping the 
distinction between gi and gf and between tu,d,e and f^^e' independently set common 



The SM as well as the MSSM limits come out right. Since an even number of SUSY particles couple at 
any vertex, we do not, however, understand how could enter without being multiplied by a second 9 
for another SUSY particle. 
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thresholds for the heavy Higgs scalars {6h) and the higgsinos (Oj^), along with a common 
threshold at rrih for the "light higgs doublet" {Oh), we find that the RGE for the coupling of 
the up-type quarks to Higgs bosons becomes, 

{Anf^^^ {3 [s'Oh + e'en] {fJlU + [c'e^ + s'Oh] (f.f]).. 



+Ac' [-Oh + Oh] (fdf])a. \ (Qkj 

4 



s 



+ 



Ikl 



kj 



{s^Oh + e'en) Tr{3fif4 + (9^ - Oh) Tr{3f]f, + fjf J 
[~9'-f {s'eh + e'en) + \~9'^f{c'eh- e'en) 
Ob [W'^\" {s'Oh + e'en) + |^"^^|' (c^^. - c'Oh)] } 



2J 



+ \9l + 8(?3' 



(36) 



Here, as with (155]) and in all the RGEs that follow, including those in the Appendix, repeated 
flavour indices k,l, . . . are summed over, including indices repeated thrice, once in a ^ and 
twice in couplings. 

It is easy to see the reduction to the MSSM, and as before, to the SM RGE for the Yukawa 
coupling matrix A„. Indeed, for Q < rriH, this RGE becomes the RGE for the coupling, 
Xu = sin /9 fu, of up-type quarks to the light Higgs doublet in (IT^ . even if higgsinos, gauginos 
or matter sfermions remain in the low energy theory. The factors of sin (3 and cos (3 in the first 
pair of curly parentheses can clearly be absorbed to turn all the f s into the corresponding 
A's. This is not, however, true in the next set of terms in square brackets where there is just 
one factor of sin/5 that combines with the (fu)j^ outside the square bracket to yield (A„).^. 
Notice, however, that the remaining couplings on this line have twiddles on them, and do not 
correspond to the quark-quark-gauge/Higgs boson couplings. Indeed, it is interesting to see 
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that we get just the right powers of sin (3 and cos (3 on the right-hand-side for all the MSSM 
couplings of matter fermions to Higgs bosons to reduce to those of the SM when all ^'s other 
than 9h vanish. Notice also that the g couplings have a non-trivial matrix structure. Since 
it is these couplings (and not the corresponding gauge boson couplings) that directly enter 
the decays of squarks and sleptons, it behooves us to keep careful track of these in a study 
of flavour-physics of sparticles. 

While it is straightforward to derive the MSSM RGEs starting from the general forms 
in (l6|)- (fTT|) . we reiterate that the reader should keep clear the distinction between the use 
of . . . as fermion field type indices in the general form from their use as flavour indices 
in the MSSM RGEs such as (1361) . Thus, in contrast to the trace in (1361) . the trace in 
refers to a sum over the fermion field types, not just flavour. This means that not all trace 
terms in fl6l)- f|TTl) lead to a trace in the RGEs of the MSSM, and further, sometimes the 
trace in the MSSM RGEs may originate in terms that are not traces in the general RGEs. 
We illustrate this with some examples. For instance, in the derivation of (!36!) . starting from 

• Since the matrix indices j, k on (UDjk and (Ua)jk denote Dirac fermion type and 
flavour, there will be a different matrix index for each quark and lepton flavour. As a 
result the term U^TrjU^^U^ + U^+U^} in (El) leads to the trace over flavours found in 
(EE]). 

• The matrix indices in (X^)jfc and (X^)jfc refer to Majorana fermions, and so do not 
carry any flavour. Tracing over these in ([6]) can, therefore, never lead to a trace in the 
MSSM RGEs. This trace term results in the term that immediately follows the trace 
term in fl36|) . 

• The situation is also different for the (Va)jk and {^a)jk in ^oi V^, flavour is 
located in the a and j indices and for in the a and k indices. When there is a sum 
over squarks or sleptons {i.e. a sum over the scalar index), we include a 6q., where i 
is a flavour index, to keep only the active sfermions at that scale. This means that we 
can keep usual matrix multiplication with the proviso that the contribution from each 
squark flavour is associated with a different 6. For example, the sum over left-handed 
squarks in 2V,X2tW, leads to -3s%^Q^^^(^'^")*(g«)f,(f«),,. 
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As a different example, consider the RGE for (g'"^) derived using (fTOi) : 



,d{i 



^ ^ dt 



Ikj 



;g"^)o- 



+3^4(g"")L(g"")''^ + 2^.-.(r«)L(r«)^^ 



^«,(g'"^)a. 



Ikl 



kl\ 



^ '2 , A 2 

3^ +%3 



(37) 

Here, the terms involving a trace over flavours may not be immediately evident. One such 
trace term, found on line four, is 



in which the trace is not explicitly written since we need to keep information about the 
position of the squark thresholds. We point out the following: 

• The terms involving the trace over flavours in fl37l) do not originate from the trace 
term in ffTOl) . The trace in fl37|) can only originate when the external quark and squark 
flavour indices on the left-hand side occur in the same or matrix on the right- 
hand side. We see upon inspection this is the case only for the term in the second line of 
fllOp . which necessarily leads to the trace over flavours in fl37l) . An analogous discussion 
applies to the derivation of the RGE for, say g"^, starting from Eq. flTTl) . Note also 
that when we are checking reduction to the MSSM, gaugino coupling matrices such as 
(g'"^) reduce to {g' x 11), so that we obtain an extra factor of three when we take a 
trace over the flavours. 
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• When inserting thetas, we must, of course, only put thetas which correspond to 
summed internal indices. To obtain, for example, the term —0i^9Q^{i^)Ji^{g,"^)ki{iu^)^jj 
which comes from 2VbW*W^ in (ITU]) , where we have picked out the left-handed 
squarks in the sum over b, we insert a 6q^, with k summed over squark flavours. Sim- 
ilarly, we insert a 6^ which comes from picking out the higgsino terms in the sum over 
Majorana fermions.^^ Finally, even though we have f^^ in our example term, there is 
no 6ix because the ur squark is the same as the one which appears on the left-hand 
side, and hence is necessarily an active squark. 

• In the 1-loop RGBs for g"^, gf^ and g^ there are no contributions from or X^, 
which is to be expected since, as mentioned before, the Higgs boson does not couple 
to the gluino. 

• Special care must be taken when considering interactions with gluinos to correctly 
evaluate the sum over colours. For example, the RGE for g"^ has some terms which 
may initially seem to lead to a trace which includes g'"«. However, the term in question 
is 

where /, k are flavour indices and c, d are colour indices. This term is clearly zero, 
since Tr{^} = 0. 

The RGEs for all the dimensionless couplings can be similarly obtained. We list these in 
the Appendix. 



V. SOLUTIONS TO THE RGES AND FLAVOUR VIOLATION IN THE MSSM 

Since our ultimate goal is to examine flavour-violation in sparticle interactions, and in 
sparticle decays in particular, we are naturally led to examine the flavour structure of various 
couplings, renormalized at the scale of the sparticle mass, i. e at scales Q typically between 
~ 100 GeV and a few TeV. 

We also sum over quark flavours, I, but there is no corresponding theta since we are always working at 
scales above mt. Therefore, all quarks arc active over the whole range of validity of our RGEs. 
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The couplings of neutral gauge bosons remain flavour- diagonal under renormalization 
group evolution even at higher loops since any inter-generation couplings would violate 
current conservation. This does not mean that their evolution is independent of Yukawa 
couplings: Yukawa couplings can, and do, enter the evolution of gauge couplings via two- 
loop contributions to the RGEs, but only as a "flavour-blind trace" . As a result of this, and 
the fact that the charged current weak interaction only couples to the left-handed fermions 
and their superpartners, inter-generation couplings in the charged current weak-interaction 
are completely determined by the Kobayashi-Maskawa 26|] (KM) matrix. 

We thus turn to an examination of the flavour structure of fermion-fermion-scalar in- 
teractions. These include the familiar Yukawa couplings of quarks with Higgs bosons, as 
well as the couplings of the quark-squark system with gluinos, charginos and neutralinos. 
Although we will conflne ourselves to the quark/squark sector, our considerations readily 
extend to the lepton sector. In this case, however, the MSSM would need to be extended to 
include effects from the flavour structure of the singlet neutrino/sneutrino mass matrix and 
interactions. 

It may appear simplest to work in a basis where the Yukawa coupling matrices in ( fT7|) are 
diagonal at one chosen scale, so that quarks are mass eigenstates, rather than in the "current 
quark basis" which we have used to write the superpotential (fT5|) and the interactions (fT6|) 
- f|T8|) as well as the RGEs. The two bases are related by the transformations, 

Vl{u)u^ = UL, Vi(rfX = cIl- (38) 
V^(mX = UR, VR{d)dj{ = dn 

where the unitary matrices \l{u) and V^r^u) [Vi^d) and VR^d)] diagonalize the Yukawa 
coupling matrix [f^] via 

Yl{u)i^VUu) = Vl{d)UVUd) = f^^ , (39) 

Of course, the matrices f^"^^ and f^'"^^ are diagonal only at one scale, which we take to be 
Q = nit. The KM matrix that enters the interactions of bosons with quarks is then 
given by 

K = Viiu)\Lid). 

The problem with choosing to work in the quark mass basis is that using different transfor- 
mations for the left-handed up and down quarks breaks the SU{2)i symmetry. As a result, 
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interactions of say the charged and neutral Higgs bosons with quarks are no longer simply 
related by an SU{2)l transformation, but involve an additional matrix (which in the case 
of the MSSM is the KM matrix). 

It is convenient to work instead in the basis where just one of up- or down- type (but not 
both) Yukawa couplings are diagonal at mj. The interactions and the RGEs that we have 
written in the last section continue to hold in this special "current quark" basis, provided we 
define the squark fields as the super-partners of these quark fields. In the following, we will 
choose the basis so that the up quark Yukawa couplings are diagonal at nit, ^-e. the up- type 
quarks are in their mass basis at this scale. This does not imply that the up-type squarks 
are also in the mass eigenstate basis. This would be the case in the exact SUSY limit, but 
the SSB squark mass matrices and the trilinear SSB scalar couplings may independently 
violate flavour in the squark sector. In writing our RGEs with independent thresholds for 
the squarks we have, however, assumed that mixing between up-type squarks is small in this 
basis so that and um are approximately also mass eigenstates, and likewise for down-type 
squarks. 

With this preamble, let us turn to an examination of the flavour violating dimensionless 
couplings of the MSSM. Let us begin with the Yukawa couplings of quarks to neutral Higgs 
bosons. We will assume that we have diagonalized the SM Yukawa coupling matrices at the 
scale Q = rrit, and examine the solutions for the various elements of the up quark Yukawa 
coupling matrix for larger values of Q. We thus need to solve the system of RGEs listed 
in the Appendix. Notice that though the RGEs for the the dimensionful SSB parameters 
are decoupled from the RGEs for dimensionless couplings, the weak scale values of the 
dimensionful parameters nonetheless enter our analysis in that they determine the various 
thresholds. This dependence is, fortunately, only logarithmic and it suffices to approximate 
the location of the thresholds as discussed in Sec. IIIII To solve the RGEs we rotate the 
Yukawa coupling matrices to the current basis using fl39l) . with matrices /j(g) chosen to 
reproduce the KM matrix. Next, we evolve these along with the values of gauge couplings 
to the scale Q = Mgut using the RGEs with all the ^'s set to unity, i.e. the MSSM RGEs. 
We then evolve back down to Q = rrit, but this time including thresholds and the couplings 

This matrix is the KM matrix for the MSSM where one of the Higgs doublets couples to up-type quarks, 
and the other to down-type quarks. In the general two Higgs doublet model, this is no longer the case. 
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with twiddles. We reset the gauge couphngs and current-basis Yukawa couphng matrices 
back to their input values, run back up to Mqut, and iterate the procedure until it converges 
with the required precision. We can then read off the couplings in any basis at the desired 
values of Q. 

Rather than show numerical results for many cases in many SUSY models that have 
been considered, we have chosen to illustrate our analysis for a simplified scenario where the 
sfermions and the electroweak gauginos are all at a mass scale ~ 600 GeV while the heavy 
Higgs bosons and gluinos have a mass ~ 2 TeV. Thus, our effective theory is supersymmetric 
for Q > 2 TeV, includes only sfermions, charginos and neutralinos together with SM particles 
for 600 GeV < Q < 2 TeV, and is the SM for Q < 600 GeV. In the intermediate range 
between 600 GeV and 2 TeV, we will have to separately examine the couplings g* and g'* 
(but not gf ) and f*^g as these split off from the usual gauge and Yukawa couplings. 



A. Quark Yukawa couplings 

The values of the running quark mass parameters determine the Yukawa coupling matrices 
in the mass basis. These can be rotated to an arbitrary current basis using fl38|) . or to the 
particular basis with diagonal up-type Yukawa couplings by choosing the four matrices 
suitably. In models where the SSB interactions are flavour-blind (at any one scale) the four 
matrices are separately unphysical, and physics is determined by just the KM matrix. In 
more general scenarios though, physical quantities will depend on more than just the KM 



matrix. In our numerical illustration, we have used the program ISAJET [27| to extract 
quark Yukawa couplings at the scale Q = rrit in their mass basis. Starting from this 



This procedure would be somewhat modified if instead of choosing the location of heavy Higgs boson 
and sparticle thresholds "by hand" as we do here, these were to be obtained using GUT scale boundary 
conditions for SSB parameters. Universal boundary conditions, such as the ones used in mSUGRA, can 
of course be used in any sfermion basis, but in all other cases we have to be specific about the choice of 
basis in which we specify these SSB boundary conditions. This will be discussed in more detail in Ref. 



2l| . As an aside, we also note that to get the solutions of the RGEs in this paper, it is only necessary to 
evolve beyond the highest threshold during the iterations to obtain convergence. We have evolved all the 
way to Mgut only because we anticipate that this will be the Q value where SSB boundary conditions 
will be specified in many realistic scenarios. 

ISAJET currently does not include flavour mixing among the quarks. It does, however, include important 
radiative corrections to the relationship between quark masses and the Yukawa couplings. We expect that 
because this flavour mixing is small, its effect on the quark masses will also not be large. In any case, we 
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FIG. 1: Evolution of the magnitudes of the elements of the up-quark Yukawa coupling matrix for 
the MSSM with thresholds, shown by arrows, clustered at 600 GeV and 2 TeV as discussed in the 
text. Above uih (= 2 TeV) we plot |(f-u)jj| whereas below rriH, where the effective theory includes 
just one scalar Higgs doublet, we plot |(AM)jj| / sin/3 which is equal to |(fM)jj| at Q = mn = 2 TeV. 



boundary condition, we have iteratively solved the system of RGEs for a toy scenario where 
the SUSY thresholds are taken to be clustered either at 600 GeV (Mi, M2, and tn^), or 
at 2 TeV {rriH, i^g), as mentioned in the previous paragraph. 

The results of this calculation for the elements of are shown in Fig. [H Specifically, 



we show 



if Q > rriH, and 



/ sin (3 (which joins continuously to the the corre- 
sponding element of at Q = ttlh) for Q < rriH- Since the KM matrix includes a complex 
phase which we take to be 60° (and since the matrices in (1551) can themselves be complex), 
these elements are, in general, complex numbers. In Fig. [1], we have plotted the absolute 
values of these elements, in the basis where the up-quark Yukawa coupling matrix is diag- 



expect to incorporate these RGEs into ISAJET, at which stage flavour mixing effects can be included in 
the iterative procedure used to extract the Yukawa couplings. Since our discussion here is meant only to 
be illustrative, it is not crucial that we use absolutely precise values for the Yukawa couplings. 
We have solved these equations in a randomly chosen current basis, i.e. one with a random choice of the 
matrices in psp , and then rotated back to the basis where up- type quark Yukawa couplings are diagonal 
at Q = rrit. Our results are independent of the initial choice of current basis as they should be, providing 
a non-trivial check on the code. 
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onal at Q = rrit, versus the energy scale Q. In our calculation, we include two-loop terms 
albeit without threshold corrections which are numerically completely negligible as even the 
two-loop corrections have only a small effect on our results. The three roughly horizontal 
lines show the scale dependence of the diagonal Yukawa couplings. This is largely governed 
by the strong interaction gauge coupling terms in the RGE, and for large Q values scales 
as a calculable power of this gauge coupling. There are kinks in these curves at 600 GeV 
and 2 TeV, but these are not visible on the scale of this figure: see, however. Fig. [51 The 
off-diagonal elements of (the scaled) |A„| all start off at zero aX Q = rrit, but rapidly increase 
to values that though small, may as in the case of {iu)23 become as large as 10~^. The most 
striking thing about the figure is the dip in all the off-diagonal elements at a value of Q in 
the few hundred GeV range, close to the expected value of squark masses in our example. 
An understanding of this feature is of more than mere academic interest, since we shall see 
that it also appears in the corresponding f-couplings, which directly enter the amplitudes 
for flavour violating squark decays. 

To understand this curious feature, let us first consider a further simplification where the 
two SUSY thresholds coalesce into a single one at Q = M. We will, of course, return to the 
case in Fig. [1] shortly. For Q > M, the RGE for f„ becomes the corresponding MSSM RGE 

^I'f diag\ 

+ (f„^-^)Tr{3(f„^-s)t(e^)} (40) 



13 2 o 2 16 2 



while for Q < M, we have the SM RGE for the Yukawa coupling matrix, 

^/XdiagN o Q 

(4vr)^^^^^ =|(Af^)(Af ^)t(Afs) - |K*(Af ^)(Af ^)tK^(Af ^) 



+ (Af s)Tr{3(Af ^)t(Af ^) + 3(Af ^)t(Af ^) + (Af ^)t(Af s)} (41) 
- (Af *^) 



The KM matrix K appears in these equations because we have written them in the basis 
with all Yukawa couplings diagonal and real at the scale Q = rrit. 

That the real or imaginary parts of the off-diagonal elements all have a zero is simple to 
understand. First, we observe that since A^ and A^ start off as real and diagonal at Q = rrit, 



18 



For the (3,3) clement, contributions from the third generation Yukawa couphng are also significant. 
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off-diagonal elements develop only because the KM matrix has off-diagonal components. 
Also, the imaginary parts of An,d develop only via the phase in the KM matrix. Next, 
notice that the "seed of flavour- violation" , i.e. the term involving the KM matrix in the 
RGEs, which is also the seed for the imaginary part of the off-diagonal elements, enters with 
opposite signs in the SM and in the MSSM RGEs. This means that if this seed term causes 
the real or imaginary part of any off-diagonal element of say the up-type Yukawa coupling 
to evolve with a negative slope in the SM, it would cause the same element to evolve with a 
positive slope in the MSSM, and vice-versa. The evolution ioi rrit < Q < M is governed by 
the SM ROE, which causes the real part of an off-diagonal element to evolve to say negative 
values. However, for Q > M, the evolution is governed by the MSSM RGE, and the real 
piece of this Yukawa coupling now begins to evolve in the opposite direction, so that it 
first becomes less negative, then passes through zero and continues to positive values. Of 
course, at Q = M, we must remember to switch from the SM couplings Xu,d,e to the MSSM 
couplings fu,d,e- The imaginary piece of the off-diagonal Yukawa couplings similarly passes 
through zero. However, we see from the sharp dips in Fig. [1] that the real and imaginary 
parts have a zero at (almost) the same value of Q, and furthermore, this location appears 
independent of the choice of the flavour indices i and j. To understand this, we need to 
analyse the equations further. 

In the following we will collectively denote the solutions of these equations by U, D and 
E, z.e U = fu for Q > M, while U = for Q < M, and analogously for D and E. For 
example, the RGE for U in the basis where both up- and down- type quark matrices are 
diagonal at Q = rrit can be written in the form, 



Here, d(t) denotes the collection of the three gauge couplings. Of course, the coefficients Ai 
and the function Q differ for the SM and the MSSM. As already mentioned, the evolution of 
the diagonal elements of U and D is dominantly governed by the last term in these equations 
which contains the gauge couplings. If we drop all other terms on the right-hand-side, the 
RGE for the diagonal Yukawa terms can be easily integrated, and we find. 



(4vr) 



dt 



+ Uij-TrjAgU^U + A^D^B + A^E^E} + U,,^(a(t)) . 



(42) 



Uii(t) = IJ ii{to)F{d{t)), no sum on i. 



(43) 
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It is important to note that the function F is independent of the quark flavour}^ 

Note now that the trace as well as the gauge coupling terms are flavour-independent. 
Assuming that the diagonal Yukawas saturate the trace (remember that we just argued that 
these diagonal terms are a function only of the gauge couplings we can combine these 

terms so that the second line on the RHS of fH21) can be written as 



u.,(t)/(t), 

again with f being independent of flavour. Next, let us look at the term involving the KM 
matrix. Noting that the diagonal elements of U and D are much larger than the off-diagonal 
ones, we can write, 

k 

where i and j are unsummed, fixed indices. In writing this, we have used the reality of the 
diagonal elements of the Yukawa matrices under renormalization group evolution (discussed 
below). From the analogue of fH5]) for D{t), we see that this term has the structure, 

M,,(to)G'(t), 

with Mij(to) = K*^Dfcfc(to)Dfcfc(to)KjfcUjj(to), and G a flavour-independent innction of just 
the gauge couplings, which themselves depend just on t. Ignoring, for the moment, the first 
term in the RGE for U (the one cubic in U) we see that we have managed to decouple the 
RGBs for Ujj which (if the first term is dropped) take the form, 

DIJ,^ = ^ = M,,{to)G{t) + /(t)U,,(t), 

where f{t) and G{t) are known functions of t. The solution to this may be written as, 

U,,(t) = M,,(to)-^G(t) + U°.(t), 



In the approximation where we retain just the last gauge couphng terms in (|40[) and (|¥T|) . these equations 
take the form, 

(4^)'^ = U,,(aig2 + a29' + a^g^). 
They can be readily integrated to give, 

u,,(^) = u„;(^o)^tl(5r)"''/'^ 

where br is the coefficient of the one-loop gauge coupling /3-function, Pridr) — j^brgf- 
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where -jj^ is the resolvent operator, and 



U°.(t)=U°(to)exp (^ffm) 



is the solution to the corresponding homogenous differential equation. For the SM evolution 
starting at Q = m^, U°j(to) vanishes (if i j), and we are left with the solution of interest, 

V,,{t) = M,,{to)n{t), (44) 

where the function 7Y(t) = j^G(t) is independent of flavour. All the flavour information 
is contained in the "boundary value", Mjj(to)- To clarify, we use the boundary conditions 
on the SM Yukawa couplings to obtain U up to the (common) SUSY threshold. The values 
of U and D at Q = 2 TeV now serve as the boundary condition for MSSM evolution. The 
central point of this analysis is that the flavour- dependence is completely captured in the 
first factor Mjj(to), where to corresponds to Q = nit for SM evolution, and to the SUSY 
threshold for MSSM evolution. 

To the extent that we can ignore the first "cubic term" in the U RGBs, we see that the 
flavour structure of the evolution of \ij is given by (l44l) . We make the following observations. 

• This structure is independent of the coefficients in the various terms in the RGEs as 
well as independent of which matrices enter into the trace. The same structure thus 
obtains for both SM and MSSM evolution. The important approximation was that 
the diagonal elements of the Yukawa matrix dominate in the basis where these are 
diagonal at the low scale. 

• Eq. ( 144|) was obtained without any assumptions about the hierarchical structure of 
the KM matrix. 

We can immediately obtain two corollaries from (l44l) . 

1. XFlfj' where the superscripts R and /, respectively denote the real and imaginary 
parts, is independent of t. 

2. —jij-j — is independent of t. 

Our program is to start at Q = rrit with the "diagonal boundary conditions" for the SM 
Yukawa couplings and evolve Xu,d,e to the scale Q = uih- We then scale the up- and 
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down-type Yukawa couplings by 1/ sin/? and 1/ cos (3, respectively, and finally continue this 
evolution to higher values of t using the MSSM RGEs. By corollary 1 above we see that 
when the real part of any off-diagonal Yukawa coupling vanishes (as it must for some value 
of Q > M), the imaginary part of this coupling also vanishes. Corollary 2 then tells us that 
the other off-diagonal elements likewise vanish for this same value of Q. This is, of course, 
exactly the qualitative feature in Fig. [1] that we started out to explain. 

Before turning to other issues, let us briefly return to the first term that we have ignored 
up to now. If we saturate the product UU^U with two diagonal elements in each of three 
terms and also remember that the diagonal elements are real, i.e. write (UU^U) ~ 
UiiUiiUij + UjjXJjjUij + UjjUjjU*j, we find that, like terms we have considered up to now, 
two of these terms do not couple Ujj with anything else, while the third term couples Ujj 
to U*j. Eq. (jUj) is thus violated only by this last term. We have indeed checked that the 
position of the zeros are not exactly coincident but have a very tiny spread of a few GeV 
not visible on the figure. We have also checked that increasing the coefficient of the cubic 
term increases this spread, while altering the other coefficients has no effect on it, in keeping 
with expectations from our analysis. 

Our analysis of Fig. [T|, up to now, has assumed a common location for the threshold for all 
non-SM particles which is, of course, not the case in the figure. It is, however, clear that the 
bulk of the change in the slope of the curve occurs from the thresholds at 600 GeV where the 
coefficient of the term involving the KM matrix changes its sign, causing the magnitudes 
of the off-diagonal couplings to decrease until they pass through zero, beyond which the 
magnitude starts increasing once again. The kink in the various curves at Q = 2 TeV marks 
the position of the second threshold. It should, therefore, be clear that the location of the 
zero is largely determined by the spectrum, and is insensitive to other details. 

Finally, we remark that the imaginary parts of the off-diagonal couplings (recall that these 

We have also numerically solved the RGEs using a fictitious KM matrix with large off-diagonal entries. 
We found the same behaviour as in the figure, in keeping with our observation that the analysis does not 
depend on the hierarchical structure of the KM matrix. 

Wc checked also that using a fictitious KM matrix with large off-diagonal entries changes the evolution 
of the off-diagonal couplings (for instance the relative size of the imaginary parts) as expected but hardly 
affects the location of the zero which is determined by the mass spectrum. As mentioned previously, the 
location of the zeros is also largely insensitive to the coefficients of the various terms in the RGEs, unless 
of course, a coefficient happens to be so large that the solution blows up before the zero can be reached. 
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have their origin in the phase in the KM matrix) may, depending on which matrix element 
we are considering, be of comparable magnitude to the corresponding real part, or much 
smaller. The relative size of the imaginary part is largely determined by the relative size of 
the imaginary part in the corresponding element of the "seed term", ~ K*DD^K-^U in the 
RGE. It is easy to see that this "seed" is real for the diagonal elements of U (or D) , which, 
therefore, continue to remain real. For the case study in the figure, the (1,3) element has 
the largest imaginary part 0(10"''), while the imaginary part of other off-diagonal elements 
is smaller by 3-6 orders of magnitude over most of the range of Q. 

We should mention that in a complete diagrammatic calculation there would be other 
non-logarithmic corrections to the couplings not included in our calculation. The corrections 
to the diagonal elements of the Yukawa coupling matrices are already included in ISAJET 

and so will be automatically included when these RGBs 



using the formulae in Ref. 
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are embedded into the ISAJET code. In contrast to a diagrammatic approach, our method 
(also the one used in ISAJET) automatically sums the potentially large logarithms that arise 
when the MSSM spectrum is significantly split. We also remark that in Fig. [T] (as well as in 
the other figures below), we have ignored any finite shifts in the coupling constants coming 
from the matching of the two effective theories (with and without the heavy particle) when 
we decouple particles at their mass scale [29|. This would potentially give jumps in the 
couplings as we cross the various thresholds which must be taken into account to achieve 
true two-loop accuracy. The RGEs that we have derived are, of course, unaffected. 



B. Higgsino and gaugino couplings to quarks and squarks 

We now turn to a study of the evolution of the couplings of gauginos and higgsinos to 
the quark-squark system. Above all thresholds, i.e for Q > rrig = mn in our numerical 
study, these are equal to the corresponding gauge and quark Yukawa couplings, but differ 
from them for the range of Q between the lower threshold at the sfermion, higgsino and 
electroweak gaugino masses and the high threshold. Of course, for Q values below this lower 
threshold, the effective theory is the SM and the couphngs g, g', f„ and cease to be 
meaningful quantitates in our simple-minded approach with thresholds being incorporated 
by step functions. Note also that because the gluino is the heaviest SUSY particle in our 
illustration, there is no gg coupling. Above all thresholds, supersymmetry relates gaugino 
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and gauge couplings so that flavour violation via gaugino interactions is not allowed. For 
Q < rriH, the connection between the gauge couplings and gaugino couplings is broken, 
and flavour-violating interactions of gauginos are no longer forbidden by any symmetry. 
Although these g*-type couplings start off proportional to unit matrices at Q = rriH, con- 
tributions to their renormalization group evolution from terms involving Higgs boson and 
higgsino couplings to quarks and squarks (see, e.g. Eq. (156!) of the Appendix), which cancel 
out if all thetas are set to unity, cause these to develop off-diagonal elements for Q < rriH- 
These off-diagonal components of g and g' (and also of the corresponding Yukawa cou- 
plings) induce ffavour- violating decays, qi — * qjZk, of squarks even if there is no exphcit 
ffavour-mixing in the SSB sector {i.e. all ffavour violation is induced via superpotential 
Yukawa couplings as is the case in many models), and hence are phenomenologically rele- 
vant. In this case, the partial widths for ffavour-violating decays will be small, so that these 
decays will likely be most important for those squarks for which flavour- conserving, tree-level 
two body decays are kinematically forbidden. Thus, the neutralino in the flavour-violating 
decay is most likely to be the LSP. Since the LSP is bino-like in many models, it would 
seem that the couplings g' would be the most relevant, but this would of course depend on 
the size of the corresponding off-diagonal component, relative to the same component of the 
other couplings. 

We begin by showing the magnitudes of the off-diagonal components of the electroweak 
gaugino- quark- squark couplings in Fig.[2]for a) the matrices g"^ and g'"^, and b) the matrices 
g*^ (there is no g"^). We obtain these from the corresponding evolution equations with the 
boundary condition that these matrices are equal to the corresponding gauge coupling times 
the unit matrix at Q = rriH- On the low Q side, we terminate these curves at the lower 
SUSY threshold below which these couplings no longer exist. Also shown for comparison 
is the (2,3) element, the one with the largest magnitude, of the up-quark Yukawa coupling 
matrix f„. Remember that we are plotting these matrix elements in the basis where the 
up-quark (but not the down-quark) Yukawa couplings are diagonal at Q = rrit. There are 
several features worth remarking about. 

1. The off-diagonal elements of the gaugino coupling matrices vary over the same broad 
range as the corresponding elements of the Yukawa coupling matrices shown in Fig. [H 
Indeed the magnitude of some of these elements considerably exceed the largest ele- 
ment of the quark Yukawa coupling matrix. It would, therefore, be dangerous to simply 
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FIG. 2: Evolution of the magnitudes of the complex off-diagonal elements of the gaugino coupling 
matrices: a) g"^ and g'"^, and b) g'^ defined in the text for the MSSM, with the positions of the 
sparticle and Higgs thresholds (again denoted by arrows) as in Fig. [TJ Also shown for comparison is 
the magnitude of the (2,3) element, the one with the largest magnitude, of the up-quark coupling 
matrix f^. The legend is in the same order (from top to bottom) as the curves. The magnitudes 
of the elements of the gaugino coupling matrices are symmetric in the flavour indices to a good 
approximation, as explained in the text. 
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disregard these when discussing sparticle flavour physics, particularly in models where 
Yukawa couplings are the sole source of flavour violation. 

2. The off-diagonal couplings of the g"^ matrix are several orders of magnitude larger 
than those of g'"«. The reason for this can be traced to the RGEs for these listed 
in the Appendix. We see from (1541) and (1561) that the evolution of the matrices g"^ 
depends on the down-type Yukawa couplings, while that of g'"-^* does not have any such 
contributions. In the basis that we are working in, these contributions are much larger 
than those from up-type Yukawa matrices (and larger than those from the gaugino 
coupling matrices which start off as unit matrices at Q = rriH), so that it is no surprise 
that the off-diagonal elements of g"^ attain larger values than those of g'"^. A similar 
remark applies to the magnitudes of the elements of g*^ shown in frame b). 

3. We also see from the figure that the magnitudes of the gaugino coupling matrices 
are symmetric under the interchange of the two indices. We have checked that this 
symmetry is not exact but (for our illustration) holds to a few parts per mille. The 
reason for the symmetry is that in the course of their evolution starting from a unit 
matrix at Q = rriH, the gaugino coupling matrices remain (approximately) Hermitian. 
This is simplest to see for the g"^ matrices, retaining only the dominant terms involving 
the down- type Higgs and higgsino coupling matrices on the right-hand-side of ( l54l) . If 
we now ignore the small difference between the Higgs and higgsino coupling matrices 
in the RGE, and keep in mind that the matrices g"^ and g'^ on the right-hand-side are 
almost the unit matrices, it is not difficult to see that the gaugino coupling matrices 
g"^ and g*^ remain Hermitian with, but only with, our approximations. A similar, 
but somewhat more involved, argument also holds for the off-diagonal elements of 
g'""^, where we have to look at all the terms in (l56l) since, in the absence of down- 
type Yukawa matrices on the right-hand-side, no single term dominates. Nonetheless, 
we have checked that if the difference between Higgs boson and higgsino coupling 
matrices on the right-hand-side can be ignored, and the gaugino coupling matrices can 
be approximated by unit matrices, the renormalization group evolution preserves the 
(approximate) hermiticity of g'"^, explaining why the magnitudes of the corresponding 
off-diagonal elements are symmetric in the flavour indices. 

Next, we turn to the off-diagonal elements of the f-matrices. We know that these will 
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FIG. 3: The real (R) and imaginary (I) part of the (1, 3) element of the up-quark Yukawa coupling 
matrix, f„, along with the corresponding elements of the matrices iu and f^^ for the MSSM with 
the spectrum as in Fig. [TJ For the solid black and red lines, we plot the elements of the matrix 
A„/sin/3 below Q = niH, as in Fig. [TJ The main figure zooms in on the low end of the range of 
Q where the Higgs boson and higgsino coupling differ from one another, while the inset shows the 
evolution all the way to Mqut- 

deviate from the elements of the corresponding Yukawa coupling matrix only below the 
threshold a.t Q = rriH, but have magnitudes similar to the corresponding Yukawa coupling 
matrix elements, whose absolute values are shown in Fig. [H It is the difference between the 
couplings of Higgs bosons and higgsinos that will be the main focus of our attention. In 
Fig. [31 we show the evolution of the real and imaginary parts of the (1,3) element of (z.) the 
Yukawa coupling matrix, f„ for Q > rriH and A„/ sin jS (which connects continuously to f^) 
for Q < niH, (ii-) the higgsino coupling matrices f^, and f^^ whose evolution is shown in 
(HSj) and (jSTI) . respectively, of the Appendix. There is no particular reason for our choice of 
the (1, 3) element (which happens to have a comparable real and imaginary piece) for the 
illustration in the figure. Here we have focussed on the lower end of Q where the and f„ 
couplings are different, while the inset shows the evolution all the way to Mqut- Several 
points may be worthy of notice. 

• For Q > rriH where the effective theory is supersymmetric, we see that the real and 
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imaginary parts of the Higgs boson and higgsino couplings separately come together 
as expected. 

• The reader can easily check that the ratio of the real and imaginary parts of the Higgs 
boson Yukawa couplings is independent of the scale as we had discussed in Sec. IV A[ 

• For Q < rriH, the higgsino couplings are split from the corresponding Higgs boson 
couplings as well as from one another by a factor of several. For instance at the scale 
of squark masses, the real (imaginary) parts of the (1,3) element of both and f"^ 
are quite different from the real (imaginary) parts of (fu)i3. It seems to us that the use 
of the evolved Higgs boson coupling in place of the corresponding higgsino coupling 
could be a poor approximation. 

• Notice that while the real and imaginary parts of (f^)i3 and (fu^)i3 come to zero at 
the same point, the position of the zero differs for the two couplings. 

Up to now we have focussed our attention on flavour off-diagonal couplings. Before closing 
our discussion, we briefly consider the effect of the thresholds on the evolution of flavour 
diagonal couplings. As an illustration, we show the evolution of the hypercharge gauge 
coupling g' and the (3,3) elements hypercharge gaugino couplings g"^ and g'"^ in Fig. HI 

The black solid line denotes the result of our complete calculation of the gauge coupling 
including both threshold and two loop effects. Also shown by green dashed and blue dotted 
lines are the corresponding results obtained at the one- and two-loop levels, but ignoring 
threshold effects. The evolution of the (3,3) element of g"^ breaks away from the evolution of 
g' at Q = mH, and is shown as the violet dot-dashed curve that terminates at Q = 600 GeV. 
The insets on the left and right show zooms of these curves at the low and high ends, 
respectively. In the inset on the left, we have also shown the evolution of the (3,3) element 
of g'"^ as the orange short-dashed line. We have checked that all the diagonal elements 
of both g"^ and g'"^ have a similar behaviour. Moreover, the (1, 1) and (2,2) elements all 
evolve essentially together with only the (3, 3) elements split from these due to top "Yukawa" 
couplings. The gauge coupling curves use the measured value of the coupling at the low 
scale as the boundary condition, while the gaugino coupling curves are obtained assuming 

■^^ In our illustration where wc have tan /3 = 4. the curves for the (1, 1) and (2, 2) elements fall in-between 
the dot-dashed violet curve for (g"^)33 and the orange-dashed curve for (g'"")33 in Fig.[31 
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FIG. 4: Evolution of the hypercharge gauge coupling at the 1-loop (green dashed line) and 2-loop 
(blue dotted line) levels (excluding threshold effects), and with the full calculation (solid black 
line). Also shown is the evolution of the (3, 3) element of g"^ (and g'"^ in the left inset). The main 
figure shows the evolution of these couplings between Q = rrit and Q = MquTj while the insets on 
the left and right zoom in on the range of Q near the TeV scale and near Mgut respectively. 

that the corresponding matrix equals the hypercharge coupling times the unit matrix at 
Q = rriH- Several points are worth noting: 

1. Since the value of the gauge coupling is fixed a,t Q = rrit in our illustration, the 
difference between one and two-loop evolution (without threshold effects) is seen at 
the high Q end where the curves are terminated at Q = Mqut defined as the point 
where the SU{2) and the scaled hypercharge couplings meet. Notice that Mgut is 
different in the two cases, as is the value of the gauge coupling at Q = Mgut- This 
relative difference is O {j^), roughly the expected magnitude of a higher loop effect. 
The full calculation of the evolution depicted by the solid black line differs most from 
the other two curves at the low end, because we start off with the evolution using 
the /3-function of the SM and, as we pass through the various thresholds, join up 
at Q = ran to the MSSM evolution. Beyond Q = rriH, the full and two-loop curves 
evolve with the same /3-function, but the latter ends up lower because threshold effects 
caused it to start off lower at Q = itih, as can be seen in the inset on the left. That it 
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ends up so close to the one-loop curve is a coincidence, but its proximity reflects the 
conventional wisdom that two loop effects in the evolution are numerically comparable 
in magnitude to the threshold corrections in the evolution at one-loop. 

2. Turning to the gaugino couplings at the low end, the most striking feature is that below 
Q = tuh, these evolve in the opposite direction to the gauge couplings. This behaviour 
can be understood if we recognize that for Q < nig the evolution of the diagonal gaugino 
couplings now depends on the much larger gluon coupling even at the one-loop level. 
As we can see from fl5^ . the terms involving cancel the terms depending on the 
QCD coupling constant Qs ii Q > rrig. However, for Q < nig the terms with 6g are 
no longer operative, and this cancellation is incomplete, causing a large change in the 
/5-function, and hence in the slope of the curve. It is striking to see that even though 
we have maintained both thresholds not far from the TeV scale, the corresponding 
gaugino and gauge couplings can develop a difference of ~ 4%. The existence of a 
difference between a gauge boson coupling and the corresponding gaugino coupling 



has been discussed in Ref. 



30| (although without any flavour structure), where it 



was suggested that its determination at an e~^e~ linear collider [3l| would give an 
idea of the splitting between sparticle masses, even if the heavy sparticles are not 
kinematically accessible. 

Finally, we turn to flavour- conserving Higgs and higgsino interactions. As an illustration, 
we show the evolution of the (3,3) elements of fu, of f^, and of f^-^ in Fig. [5l Just as in the 
previous flgure, we show results for the complete calculation at the two-loop level, including 
threshold effects as well as differences between the couplings of Higgs/gauge bosons and of 
higgsinos/gauginos by the solid black line. We also show the results that we obtain with 
just MSSM evolution all the way to Q = rrit using one-loop RGEs (the green dashed curve) 
and two- loop RGEs (the blue dotted curve). In all three cases we start with the same value 
for the Yukawa coupling a.t Q = rrit. Also shown by the violet dot-dashed and orange short- 
dashed lines is the evolution of the (3,3) element of the matrices and f^^, respectively. 
The two insets are similar to those in the previous flgure. We note the following: 

1. We see that the complete calculation of the Yukawa coupling leads to a large difference 
from the two-loop calculation without thresholds over most of the range of Q. This is 
not new and is largely due to the difference between the evolution of Yukawa couplings 
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FIG. 5: Evolution of the (3,3) element of the Yukawa coupling, £„, at the 1-loop (green dashed 
line) and 2-loop (blue dotted line) levels (excluding threshold effects), and with the full calculation 
(solid black line). Also shown is the evolution of the (3, 3) element of (and f"^ in the left inset). 
The main figure shows the evolution of these couplings between Q = rrit and Q = MquT) while 
the insets on the left and right zoom in on the range of Q near the TeV scale and near Mqut 
respectively. As in Fig. [H we have plotted A^/ sin/3 below Q = mu- 

in the SM and in the effective theories that interpolate between the SM and the 
MSSM.^^ Below the kink at Q = 600 GeV in the solid black curve the evolution of the 
Yukawa coupling is as in the SM, while above Q = rriH it is as in the MSSM. Notice 
that in this case, threshold effects are considerably larger than the difference between 
one- and two-loop evolution. 

2. It can be clearly seen from the left-inset that the higgsino couplings evolve quite 
differently from the corresponding Higgs boson couplings once Q is below rriH- Once 
again, this difference is largely due to incomplete cancellations in terms involving 

Notice in the left inset that the solid black line "curves" significantly between Q = 600 GeV and Q = mt- 
This highlights the advantage of this approach (also used in the event generator ISAJET) which "sums 
logs of the ratio of the low and high threshold by solving the RGEs", with that sometimes used in the 
literature where MSSM evolution is used all the way down to Q = mt and then corrected for "via a single 
step" evolution to take into account the difference between the running in the MSSM and in the SM. 



50 



"strong" interaction couplings of gluons and gluinos, as may be seen from (HHjl and 
fl5T]) of the Appendix. 

VI. CONCLUDING REMARKS 

RGEs provide the bridge that allows us to extract predictions of theories with simple 
physical principles valid at very high energy scales, many orders of magnitude larger than 
energies accessible in experiments. Because of renormalization effects, these same simple 
principles lead to complex predictions at accessible energies. Since supersymmetric theories 
allow sensible extrapolation to high energy, RGEs have played a central role in the analysis 
of many supersymmetric models, generally assumed to reduce to the MSSM (possibly aug- 
mented by right-handed neutrino superfields) in the range between the weak and GUT or 
Planck scales. 

In this, the first of a series of two papers, we have re-examined the threshold corrections 
to the RGEs for the dimensionless couplings of the MSSM, incorporating also the effects 
from flavour mixing of quarks and squarks. Above the scale of all new particle thresholds, 
the effective theory is the MSSM, with just three gauge couplings and three different Yukawa 
coupling matrices (that specify the interactions of matter fermions with Higgs bosons) being 
the independent dimensionless couplings of the theory. All other dimensionless couplings, 
for instance those of gauginos or higgsinos to the fermion-sfermion system, or quartic scalar 
couplings, are related to these by supersymmetry. These relations are, however, no longer 
valid once supersymmetry is broken, so that then the couplings to gauginos and higgsinos 
will renormalize differently from gauge and Yukawa couplings, respectively. In a consistent 
treatment of threshold corrections, the RGEs for these couplings will therefore differ from 
the RGEs for the gauge and quark Yukawa couplings that are available in the literature. 

We have adapted the RGEs for a general {i.e. non-supersymmetric) field theory Q, [3] 
and rewritten them in four-component spinor notation that we use when we obtain the 
RGEs of the MSSM. The details of our procedure may be found in Sec. |TT] and llllt and 
the application to the MSSM in Sec. IIVI The complete set of RGEs for the dimensionless 
couplings between SM particles and their superpartners (but not for the quartic couplings 
of scalars) is given in the Appendix. These quartic couplings are less important from a 
phenomenological perspective, though some of these do enter the squark and slepton mass 
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matrices. The important thing for the present discussion is these do not enter into the RGEs 
for the fermion-fermion-scalar couphngs, which can then be evolved independently from the 
quartic couplings. 

We have presented some sample numerical results for the evolution of various couplings in 
Sec. El The analysis here is meant only to give a sampling of effects that, to our knowledge, 
have not been previously pointed out in the literature. Since flavour physics in the sparticle 
sector has been the main motivation for our analysis, we point out that when threshold effects 
are included, gauginos are not different from higgsinos in that neutral gauginos also develop 
flavour- changing couplings to quarks and squarks. As seen from Fig. O while these flavour- 
violating gaugino couplings vary over several orders of magnitude, they can be comparable, 
or even larger, than the corresponding couplings to Higgs bosons. An illustration of flavour- 
violating higgsino couplings is shown in Fig. [31 The scale-dependence of these is quite 
different from that of the usual Yukawa couplings of quarks to Higgs bosons. We stress that 
the induced flavour violating couplings of gauginos and higgsinos, evaluated at the scale of 
any particular squark mass, will contribute to the amplitude for the flavour-violating decay 
of this squark. Such a contribution is distinct from the usually evaluated contribution from 
the induced flavour-violation in the squark mass matrix. Both types of contributions need 
to be included. Finally, we point out one last striking feature of our analysis. We see from 
Fig. m and Fig. [51 that the inclusion of SUSY threshold corrections leads to a difference of 
a few percent between the flavour- diagonal electroweak couplings of gauginos and gauge 
bosons, and of higgsinos and Higgs bosons, traced to SUSY QCD contributions in even the 
1-loop RGEs. This could have a noticeable effect on the evaluation of sparticle masses as 
well as of thermal relic densities of neutralino dark matter created in the Big Bang. 

In summary, we have presented the RGEs of the dimensionless couplings of the MSSM 
in this paper. In a follow-up paper under preparation 2]|, we will present the RGEs for 
the dimensionful parameters, again including threshold and flavour-mixing effects. This 
complete set will facilitate the examination of flavour phenomenology in SUSY models with 
arbitrary ansatze for flavour in the SSB sector. 



52 



Acknowledgment s 



We are grateful to H. Baer, D. Castafio, A. Dedes, S. Martin, K. Melnikov, A. Mustafayev 
and M. Vaughn for clarifying comments and communications. We thank D. Castafio for 
sending us the unpublished erratum to Ref. 1^, and H. Baer and A. Mustafayev for their 
comments on the manuscript. This research was supported in part by a grant from the 
United States Department of Energy. 

Appendix 




53 



(47r) 



dt 



+4^2 [_^, +^^] (f„ft),,} (f,),^. 



+ 



C(fd)jfc 

c 



■>s J Ij 



Ikl 



Qfclg?)iA;(g< 



2 
3 

+ c(fd)ii 

c 

+ 5 



{Oh - Oh) Tr{3ftf4 + (c^^, + s^^^) Tr{3f]f, + ftfe}" 



- c(f.) 



^^7" + ^^72 + 8^?' 



(46) 



(47r 



,2'^(cfe)i 



---c {cHn + s^e^,] (f,ftf,),, 

+ c(f,),, [^^,%,(fe^)L(fe^)u + e^^.,(r«)fcz(r«)5 



c 

+ 4 



'x.'L\T f~iL\* 



2^;,^e-,(f|«)a-(f:-)L + 3^^^^^^(g^)4(g^)I, + ^^^^^(g' 



+ ce-,h, [-3^vF(^'^)*(g")f. + ^i^(^''^)*(g'^)a (fe^)/^. 

-2ce^^^e,,(^"^^)*(f:«),,(r^)j,. 



+ c(fe),, (^, - Tr{3ftf4 + {c^Or, + s^^^) Tr|3f]f, + fjf, 



c 

-t 

2 



C(fe)i 



1 5 /o 9 o 

+ 4^' 



(47) 



54 



(4vr) 



dt 



ik 



9 



+ \h + c'^^] (f.^).. {3^^ \t-t + 



^e^(g'3)*,(g«)r, + ^^^(g'«)*,(g'«)f, + ^%(g?):.(g?; 



(f.^)..(fn^)L + (f?).;^(f?)L (fn^) 



25 9 



(48) 



)lk 



+3^Q,(f?)L(f?)^. + ^L,(fe")L(f. 



+ \h [c'Oh + 3^^^] (f,^)., {3^,^ l^^'^f + 



16 



T 



-^i^^j,(g'«)*.(f^);cKg"^^)5 - y^5^j,(gs^)*.(ff );cKg^;., 



(49) 



55 



dt 



'B^e;V& Jkl\o / ij ■ - n- i^i 



-2Me-,(r")*.(f|«).Kg"^)S 



ik\ 



)kl\ i^e)lj 



+ %^L(fe^)a.(fi^)L(f^)/.-(f. 



e JiJ 



15 



.'2 , 9 2 



4^2 



+ 4% [^'^/^ + c'9h] (f^)., {3^^ + 9^ \f^-\' 



+ 



+ I^^B^Q.(g"')r^(g"')M + 3%^Q,(gs^)f.(g?)*.. 



^^B^Q..(g'^)f.(f.^)H(g"«)r, - y%^Q,(g?)f.(f«^);^Kgr)r, 



'5 



9 



56 



(4vr) 



dt 



^3— 



e Jkl 



+ ]0h [c'Oh + s'9h] (f^)., {39^ \~g'^\" + 9s \~9"' 



1 

+ 2 



16 



3%(gf^)r.(g^)r, 



(47r 



+ [c'dH + s'^h] (f|-)., l/'^r + \9^'\' 



1 

-i 

2 



;ifcv& )ki 



e 



If + ^sl 



57 



(4vr) 



4^. 



+ 7T^4(g"«)/.( 



lk\ 



ikl 



2 ^' 



1 1 

)L + fBis'%k{g"')l + y%(g?).;.(g?)L 



/2 , 2 , /I 2 



(54) 



(4vr: 



2C?(g 



dt 



^Q,(g'«),,(g'«)L + ^,, 



;g'")..(g'^)L 



lk\ 



+ ^%^B(g'^).. { [^'^^ + c'Oh] |/-|' + [c^^, + s'Oh] \~g 



2 

+ ^(g"').fc 



2 



3^v^(g^)..(g^)L + ^B(g'")..(g'")L 

3 



+^A(fe):.(fe)..(g'")..-(r 



Q 

/2 , 2 

4^ +4^2 



(55) 



58 



(47r 



dt 



(47r 



dt 



kl 



kl 



■;i.'UR 



)0- 



-^^(g'"-)L(g'"«)., + -%(gr)L(gr)z, 



(56) 



+ 



9 



fQ\T (fQ\* 



'^h^qS^d )iki^d )kl 



1 

-( 

2 



+ ^^i^(g"«) 



^%,(g'^)L(g'^)/. + ^L.(r)L(r)^.^ + ^^..(g'"«)L(g"") 



2 



2^^%,(f.^)f.(g'^)H(f^)?, + 2^;-^4(r^).,(f,^-)rz(f.'«);, 



+ ^4(g"")«'^ 



g^B(r«)L(r-).,+ 3^p(g.'«)L(gf^)/. 



(57) 



59 



(47r 



dt 



kl 



+ loMs'-^'h^ { [s'Oh + e'en] + [c'e, + s'e^] } 



(58) 



2 
1 

+ 



~Jhu [(,2^^ ^ ^2^^^ Tr{3ftfJ + (0, - ^- ^ 



+ 55- 
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^ ^ dt 4 ^ 
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e J ml 



+ eg"''' 
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dt 
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2 - n- VK(g*^)*i 
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iky& )ki\ IS -jij 
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1 1 fi 
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